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ABSTRACT 

We perform a systematic study of generic accidental Higgs-family and CP symmetries that 
could occur in the two-Higgs-doublet-model potential, based on a Majorana scalar-field for- 
malism which realizes a subgroup of GL(8, C). We derive the general conditions of convexity 
and stability of the scalar potential and present analytical solutions for two non-zero neutral 
vacuum expectation values of the Higgs doublets for a typical set of six symmetries, in terms 
of the gauge-invariant parameters of the theory. By means of a homotopy-group analysis, 
we identify the topological defects associated with the spontaneous symmetry breaking of 
each symmetry, as well as the massless Goldstone bosons emerging from the breaking of the 
continuous symmetries. We find the existence of domain walls from the breaking of Z2, CPl 
and CP2 discrete symmetries, vortices in models with broken U(l)pQ and CPS symmetries 
and a global monopole in the SO(3)jjp-broken model. The spatial profile of the topological 
defect solutions is studied in detail, as functions of the potential parameters of the two-Higgs 
doublet model. The application of our Majorana scalar-field formalism in studying more 
general scalar potentials that are not constrained by the U(1)y hypercharge symmetry is dis- 
cussed. In particular, the same formalism may be used to properly identify seven additional 
symmetries that may take place in a U(l)Y-invariant scalar potential. 



PACS numbers: 11.30.Er, ll.30.Ly, 12.60.Fr 



Contents 

1 Introduction 

2 Two Higgs Doublet Model Potential 

2.1 Convexity and Stability Conditions 

2.2 The Majorana Formalism 

2.3 The Vacuum Manifold 

3 Neutral Vacuum Solutions of the HF Symmetries 

3.1 Z2 Symmetry 

3.1.1 Neutral Vacuum Solutions from a Singular Matrix N 

3.1.2 Z2 Topology 

3.2 U(l)pQ Symmetry 

3.2.1 Neutral Vacuum Solutions from a Singular Matrix N 

3.2.2 U(l)pQ Topology 

3.3 S0(3)hf Symmetry 

3.3.1 Neutral Vacuum Solutions from a Singular Matrix N 

3.3.2 S0(3)hf Topology 

4 Neutral Vacuum Solutions of the CP Symmetries 

4.1 CPl Symmetry 

4.1.1 Neutral Vacuum Solutions from a Singular Matrix N 

4.1.2 CPl Topology 

4.2 CP2 Symmetry 

4.2.1 Neutral Vacuum Solutions from a Singular Matrix N 

4.2.2 CP2 Topology 

4.3 CP3 Symmetry 

4.3.1 Neutral Vacuum Solutions from a Singular Matrix N 

4.3.2 CP3 Topology 



5 Topological Defects in the 2HDM [49 

5.1 Domain Walls [50 

5.1.1 Z2 Domain Walls M 

5.1.2 CPl Domain Walls M 

5.1.3 CP2 Domain Walls M 

5.2 Vortices l60 

5.2.1 U(l)pQ Vortices leO 

5.2.2 CPS Vortices m 

5.3 Global Monopoles |65 

5.3.1 S0(3)hf Global Monopoles l65 

6 The U(l)Y-Violating 2HDM m 

7 Conclusions [72 
Appendices [74 
A cr'' Matrix Identities [74 

B The Form of and the Transformation Matrices [75 

B.l The U(1)y Constraint on S^^ [tS 

B.2 The Majorana Constraint on |76 

B.3 The Majorana Constraint on GL(8, C) [77 

C Trace and Determinant Relations for N^,^ and L^,, [78 

D Inverting the Transformation Matrix Relations |82 



1 Introduction 



The standard theory of electroweak interactions, the Standard Model (SM) [H [21 [3], is a 
renormahzable theory with a minimal particle content which realizes the Higgs mechanism 
in El m [7] to account for the origin of mass of the charged fermions and the and Z 
bosons. The SM describes the experimental data collected over the years at the LEP collider, 
Tevatron and in a number of low-energy experiments with remarkable success [8j . In spite of 
its conspicuous success, however, several key questions remain unanswered within the SM, 
such as the stability of the gauge boson masses under quantum corrections, the possible 
unification of the strong with the electroweak forces, the Dark Matter problem and the 
existence of new sources of CP violation to account for the observed baryon asymmetry in 
the Universe. 

Supersymmetric theories softly broken at the TeV scale provide a natural framework to 
successfully address all the above problems (for a recent review, see [9]). In particular, the 
Minimal Supersymmetric extension of the Standard Model (MSSM) requires the existence 
of one more Higgs doublet 02 in addition to the SM Higgs doublet 0i, so as to maintain the 
holomorphicity of the superpotential and ensure the cancellation of the chiral anomalies. In 
the MSSM, CP-even [TOl ttH [12] and CP-odd |l3l [H [I5] radiative corrections to the scalar 
potential can be very significant, giving rise to an effective CP- violating potential [T6 | [T7 t \T8 \ 
[19] which acquires the form of the Two- Higgs Doublet Model (2HDM) [20] 0. 

Recently, a classification of all the possible accidental symmetries that could occur 
in a 2HDM potential has been attempted [22|, |23l [2ll [25]. Such a partial classification 
was motivated by the use of a gauge-invariant bilinear scalar-field formalism based on the 
SL(2,C) group [2S1[2ZI[2S|, or its SU(2) subgroup [301 [22] §• The latter subgroup emerges 
as a reparameterization group of the 2HDM potential [32| in the restricted two Higgs-doublet- 
field basis 0i^2 upon canonical renormalization of possible loop-induced Higgs-mixing ki- 
netic terms [23! • detail, the 2HDM potential may exhibit accidental symmetries, for given 
choices of its theoretical parameters, and following the terminology in [231 [25], there exist 
two classes of symmetries. The first class of symmetries involve the transformation of the 
two Higgs doublets 0i 2, but not their complex conjugates 0^ 35 are called Higgs Family 
(HF) symmetries. The second class linearly maps the fields 0i^2 into their CP-conjugates 0^ 2 

^Historically, the bilinear mass operator, (rn\24>\4'2 + H.c), was missing in the original article by 
T.D. Lee |20) . However, it is worth mentioning that this dimension-two operator plays an important role in 
the renormalization of the general 2HDM potential [H], including the renormalization of possible CP-odd 
tadpole graphs [T?l . 

^Note that the largest possible symmetry group of the 2HDM is 0(8) [Slj, giving rise to a large number of 
symmetry breaking patterns, beyond the restricted set considered so far which realize 0(3) and its maximal 
subgroups. 

■^As we will see, however, the maximal reparameterization group of the 2HDM potential is GL(8, M), which 
acts on the 8 real scalar fields contained in the two Higgs doublets (t>i,2 and includes gauge transformations. 
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and are therefore termed CP symmetries. 

Three physically interesting HF symmetries of the 2HDM that have been discussed 
extensively in the literature are: the Z2 discrete symmetry [31], the Peccei-Quinn symmetry 
U(l)pQ [35] and the HF symmetry S0(3)hf [311 [221 El |25] which involves an SU(2)hf/Z2 
rotation of the Higgs doublets 0i^2- Likewise, three typical CP symmetries of the 2HDM 
that have received much attention are: the CPl symmetry, which realizes the canonical CP 
transformation 0i(2) — > 0*^2) [201 ISH ISH] , the CP2 symmetry, where 0i(2) — > (— )02{i) [33; 
the CP3 symmetry, which combines CPl with an SO(2)jjp/Z2 transformation of the fields 

01,2 [221 [231 [2a Eg. 

In this paper, we introduce a Majorana scalar-field basis where both the HF and CP 
symmetries can be realized by acting on the same representation of Higgs fields. To this 
end, we extend the aforementioned gauge-invariant bilinear formalism to the larger complex 
linear group GL(8, C), which is then reduced by a Majorana constraint and gauge invariance. 
Specifically, GL(8, C) is the reparameterization group acting on the 8-dimensional complex 
field multiplet $ that contains the two Higgs doublets 0i,2 and their hypercharge conju- 
gates ia'^(f)l 2 cLS components, where cr^ is the second Pauli matrix. The multiplet $ satisfies 
the Majorana constraint $ = C$* which, together with the constraint of SU(2)l (g) U(1)y 
gauge invariance, reduces GL(8,C) into two subgroups isomorphic to GL(4, M), where C is 
a charge-conjugation matrix defined in Section 2. The first subgroup is related to the HF 
transformations and the second one to the generalized CP transformations on the Majorana 
field multiplet Therefore, we refer to the above description as the Majorana scalar-field 
formalism, or in short, the Majorana formalism. 

As we will explicitly demonstrate in Section [HI the GL(8, C) Majorana formalism has 
the analytical advantage that scalar potentials being only constrained by the SU(2)l gauge 
group, but not by U(1)y, can be described in a similar quadratic form as in the usual 
SU(2)l ® U(1)y gauge-invariant 2HDM. In particular, the same formalism can be used to 
identify symmetries of U(l)Y-invariant 2HDM potentials that are larger than 0(3) in the 
bilinear field space, such as 0(8) and 0(4) (g) 0(4) in the real field space [3T]. As we will see 
in Section [21 these latter symmetries fail to be captured by the restricted framework of the 
SL(2,C) bilinear approach adopted in the recent literature. 

In this article, we also derive the complete set of algebraic conditions for the convexity 
of the general CP-violating 2HDM potential and its boundedness from below, by applying 
Sylvester's criterion (see, e.g. [M])- These algebraic conditions extend previous partial results 
obtained in the literature for particular forms of the 2HDM potential [311 [321 [ID] ci-nd may have 
a geometric interpretation in terms of conical sections as presented in [27]. Following [23121], 
we employ the Lagrange multiplier method to analytically calculate all non-zero neutral 
vacuum expectation value (VEV) solutions for the Higgs doublets 0i,2, associated with the six 
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generic HF and CP symmetries mentioned above. The non-zero VEV solutions are expressed 
entirely in terms of the gauge-invariant parameters of the theory, thereby obtaining the 
analytical dependence of possible non-trivial topological features in the vacuum manifold. 
As a cross-check, we verify that our solutions satisfy the minimization conditions derived by 
more traditional methods as explicitly given, for example, in [TB] . 

In order to get a topologically stable solution in the 2HDM, both the VEVs of the two 
Higgs doublets should be non-zero, such that the topological configuration cannot be removed 
away by SU(2)l ® U(1)y gauge transformations. We use an homotopy-group analysis to 
determine the nature of the topological defects associated with the spontaneous symmetry 
breaking of each symmetry. More explicitly, topological defects, such as domain walls, strings 
or vortices and monopoles, are created, when a symmetry group G of the Lagrangian, which 
may be either local, global or discrete, breaks down into a subgroup H, in a way such that the 
vacuum manifold Ai = G/H is not trivial. Knowing the topological properties of the vacuum 
manifold M. under its homotopy groups, n„(A^), determines the nature of the topological 
defects [311112]. Thus, domain walls arise for no(A^) 7^ I, strings or vortices for Ili{J\4) 7^ I, 
monopoles if 112 (A^) 7^ I and textures if n„>2(A^) 7^ I [HI 133], where I is the identity 
element. After having identified the precise nature of the topological solution, we then study 
quantitatively their spatial profile, as a function of the potential parameters. The results of 
our analysis may be used in future studies to derive cosmological constraints on the 2HDM, 
or on inflationary models with related SU(2) group structure [331 133 I36l 133 I3H|- 

The layout of the paper is as follows. Section 2 briefly reviews basic aspects of a general 
tree-level 2HDM potential, on which we derive the sufficient and necessary conditions for its 
convexity and boundedness from below. In the same section, we introduce the Majorana 
scalar-field formalism for describing the 2HDM potential, as well as possible extended scalar 
potentials that are not constrained by the U(1)y hypercharge group. In addition, we present 
the group structure of the vacuum manifolds for a set of six generic HF and CP symmetries. In 
Section 3, we calculate the neutral vacuum solutions for the three HF symmetries, Z2, U(1)pq 
and S0(3)hf, and identify their topological properties. Correspondingly, Section 4 discusses 
the neutral vacuum solutions and their topology, for the CP symmetries: CPl, CP2 and 
CP3. In Section 5, we perform a quantitative analysis of all the topological solutions found 
above, in terms of the fundamental parameters of the theory. We present several key features 
of the topological defects, including their spatial profile and energy density. In Section [6l 
we show how the Majorana scalar-field formalism can be extended to study U(1)y- violating 
2HDM potentials. We also show how the same extended version of the formalism can be 
used to identify further accidental symmetries that could take place in a U(l)Y-invariant 
2HDM potential. Finally, Section 7 contains our conclusions. Some technical details of our 
analytical calculations are presented in Appendices A-D. 
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2 Two Higgs Doublet Model Potential 



In this section we first review the 2HDM potential in the bihnear field formalism [251 123 122] • 
We then derive the conditions for convexity and stability of the general 2HDM potential, and 
briefly explain the Lagrange multiplier method for flnding the neutral VEV solutions for the 
two Higgs doublets. We then proceed by introducing our Majorana scalar-field formalism and 
present the group structure of the six generic symmetries that may occur in the 2HDM po- 
tential. Finally, we discuss the general group-theoretical properties of the vacuum manifold, 
which enable us to identify the exact nature of the topological defects in the 2HDM. 

Let us start our discussion by writing down the tree-level structure of the general 2HDM 
potential V: 



V 



— m 



2 



^I02f + f( 



(2.1) 



It is easy to see that the 2HDM potential V contains 3 mass parameters fif, ^\ and 
and 7 quartic couplings Ai,2,...,7- For the potential V to be Hermitian, the parameters ^\ 2 
and Ai^2,...,4 are constrained to be real, whereas ^^'^ -^5,6,7 ^ire in general complex. In 
order to evaluate the VEVs of the Higgs doublets 0i and 02, we have to calculate first the 
extremization conditions by solving the two coupled cubic equations 



dV 



+ 



dV 

d4 



+ 



-fil + 2Ai(0l0i) + X^icPlh) + A6(0l02) + X;{44>^) 

-ml^ + A4(050l) + A5(0l02) + A6(0l0l) + Xj{4(p2) 

-fil + 2X2{4<p2) + A3(0l0i) + A7(0l02) + x;{4<p^) 

-mil + A4(0l02) + A*(0^0i) + A*(0l0i) + A;(0502) 



(2.2a) 



. 



(2.2b) 



Finding analytical solutions to the above coupled cubic equations for the VEVs of 01^2, in 
terms of the gauge-invariant potential parameters, is a formidable task within the 2HDM. 
This problem is usually avoided in the literature, by assuming that the VEVs of 0i 2 are the 
input parameters, for a given set of quartic couplings, whereas the potential mass parameters 
are derived from these (see, e.g. [H]). Nevertheless, it would be highly preferable, particularly 
in the study of topological defects, to devise a method, in which the VEVs of 0i,2 can be 
analytically expressed, in terms of the gauge-invariant mass terms and quartic couplings of 
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the 2HDM potential. 

An analytical method which can address this problem is the bilinear scalar-field for- 
malism introduced in f2E\ f27\ . According to this formalism, the 2HDM potential V given 
in ( 12. ip can now be expressed in full by the 4-dimensional vector 



v 



^102 + 0^01 
')\(f)2 - 0201 



- 0202 



(2.3) 



where = (0i , 02)""" and (with = 0, 1, 2, 3) denote the two-by-two identity and the 
three Pauli matrices: 



a 



1 
1 



a 



1 

1 



cr 




a 




(2.4) 



It is obvious that the scalar-field multiplet spans an SL(2, C) group space similar to the 
spinorial Weyl space. Hence, the vector R'^ becomes a proper 4-vector in the Minkowski 
space, described by the fiat metric r/^,^ = diag(l, —1, —1, —1). In terms of the 4-vector R'^, 
the 2HDM potential reads: 



V 



where M„ and L^j, are given by 



IM^R^ + ^L^.R'^R'^ + Vo 



(2.5) 



M,, 



( /i? + , 2Re(m?2) , -2Im(m?2) , - /^2 ) . 

f Ai + A2 + A3 Re(A6 + A7) -Im(A6 + A7) Ai - A2 \ 

Re(A6 + A7) A4 + Re(A5) -Im(A5) Re(A6 - A7) 

-Im(A6 + A7) -Im(A5) A4 - Re(A5) -Im(A6 - A7) 

y Ai - A2 Re(A6 - A7) -Im(A6 - A7) Ai + A2 - A3 J 



(2.6a) 



(2.6b) 



Notice that we have added a constant term Vq to the scalar potential V in (12. 5p . which is 
adjusted such that the minimum of the potential Vmin is set to zero, thereby accounting for 
the vanishing small cosmological constant. 
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2.1 Convexity and Stability Conditions 



An obvious advantage of the bilinear scalar-field formalism is that the 2HDM scalar potential 
V in f l2.ip has been reduced from a fourth order polynomial in 0i 2 to a polynomial of second 
degree in R^, as given in fl2.5p . We can now calculate the neutral vacuum solutions of the 
potential V(R'^), which amounts to finding the local extrema of V(R^), for which R^ is a 
null vector, i.e. = R'^R^ = 4(0|0i)(0202) — 4:{4>\4'2){4'l4'i) = 0. To enforce the null norm 
restriction on R'^, we introduce the Lagrange multiplier ( and modify the potential V of (12. 5p 
to 

Vc = - ^M^R^ + ^N^.R'^R^ + Vo , (2.7) 

with N^,^ = L^i, — (rj^^. More explicitly, the modified quartic-coupling matrix N^,^ is given 
by 

/ Ai + A2 + A3 - C Re(A6 + A7) -lm(A6 + A7) Ai - A2 \ 

Re(A6 + A7) A4 + Re(A5) + C -Im(A5) Re(A6 - A7) 

-lm(A6 + A7) -Im(A5) A4 - Re(A5) + C -Im(A6 - A7) 

y Ai - A2 Re(A6 - A7) -Im(A6 - A7) Ai + A2 - A3 + C / 

(2.8) 

Consequently, within the bilinear scalar-field formalism, the extremization conditions for the 
neutral vacuum solutions of the 2HDM potential are given by dV(^/dIi^ = and dV(^/d( = 0, 
or equivalent ly by 

= N^,R\ (2.9a) 
R^R'^ = . (2.9b) 

For an extremal point to be a local minimum, we require that the Hessian H derived from the 
scalar potential V(0i^2) be positive definite. The Hessian H is, in general, an 8 x 8-dimensional 
matrix obtained by double differentiation with respect to all 8 scalar fields contained in the 
two Higgs doublets 0i,2, evaluated at the neutral VEVs f°2 of 4>i,2 and their possible relative 
phase ^ (for exact notation, see Section [273]) . However, for the given HF and CP symmetries, 
it is sufficient to examine the positivity of H derived in the restricted 3-dimensional space 
of and ^. Having identified all local minima, we then compare the values of the 2HDM 
potential V at these minima. The lowest value obtained for V singles out the global minimum, 
provided V itself is bounded from below. It is therefore important to derive the constraints 
on the theoretical parameters for having a scalar potential which is convex and therefore 
bounded from below. To ensure this, we require that the matrix L^j^ be positive definite [27] . 
The latter can be enforced by applying Sylvester's criterion which yields the following general 
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restrictions: 



Ai + A2 + A3 > , 

(Ai + A2 + A3)(A4 + i?5)-(i?6 + i?7)' > 0, 



(2.10a) 
(2.10b) 



(Ai + A2 + \3)i\l - lAsH - A4 [(i?6 + Rrf + {h + hf] 
-2/5 (i?6 + Ri) ih + h) + [{R& + R-i? - ih + /t)'] > . 



(2.10c) 



In the above, we used the shorthand notation: Rk = Re(Afc) and Ik = Im(Afc). In addi- 
tion to fl2.10ap - fl2.10cl) . we require that the determinant of L^,^, which is given analytically 
in f lC.Gip . be positive as well, i.e. det [L^j,] > 0. 

We may now observe that if R'^R^ > 0, this would imply that C = 0- This is since the 
2HDM potential should not modify by the addition of the Lagrange multiplier (, i.e. = V. 
Hence, possible solutions with C = usually signify a charged-breaking vacuum for the six 
HF/CP symmetries considered here and they are therefore rejected in our analysis. As a 
consequence, there are two distinct sets of U(l)gjj^-preserving minima that could occur in the 
2HDM, depending on whether det[N^;^] vanishes or not. If N^,^ is not singular, i.e. det[N^y] 7^ 
0, the vector R'^ can be obtained by simply inverting fl2.9ap . i.e. 



As we will see in Sections 3 and 4, the neutral vacuum solutions for the generic HF and CP 
symmetries under study (with exception of the CPl symmetry) imply that at least one of 
the VEVs of 01,2 is zero, when det[N^,^] 7^ 0. Such vacuum solutions are uninteresting, since 
they do not lead to stable topological defects. 

The second set of neutral vacua occurs, when the modified quartic-coupling matrix N^j, 
is singular, i.e. when det[N^y] = 0. In this case, the Lagrange multiplier ( takes on a specific 
value which leads to a singular matrix N^,^. If this happens, the undetermined component 
of R^ is calculated by requiring that the neutral vacuum condition R^R^ = is met. In this 
second class of solutions, both the VEVs of the Higgs doublets can be non-zero, leading to 
the interesting topological solutions which we study. 

For each of the neutral vacuum solutions we obtain by the Lagrange multiplier and 
Hessian methods outlined above, we cross-check that they also satisfy the convexity and 
the conventional extremization conditions fl2.2ap and f]2.2bp . In this way, we ensure that a 
stable and global neutral vacuum is found for the 2HDM potential. Since the matrix N^j, 



(2.11) 



and the Lagrange multiplier must guarantee that R'^R^ 



0, i.e. 




(2.12) 
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plays an instrumental role in our analysis, Appendix O contains analytical expressions for its 
determinant, as well as solutions for the Lagrange multiplier ( that give rise to a vanishing 
determinant, i.e. det[N^j,] = 0. 



2.2 The Majorana Formalism 

It would be interesting to introduce a formalism where both the HF and CP symmetries can 
be realized by acting on the same representation of scalar fields. For this purpose, we extend 
the gauge- invariant bilinear formalism based on the SL(2, C) group to the larger complex 
linear group GL(8, C) (see also [39] for a related discussion). Specifically, this latter group is 
acting on the 8-dimensional complex field multiplet 



\ 



(2.13) 



Notice that under a SU(2)l gauge transformation Ul, all doublet components of the multiplet 
$ transform in the same way, i.e. $ — )■ Ul with 



Ut 



exp 



16' 



a 



0-° (g) cr° ® exp 



«eVV2 



(2.14) 



where the summation convention over the repeated group indices i = 1, 2, 3 is assumed, with 
^1,2, 3^-2 being the generators of the SU(2)l gauge group and 9^'^''^ G [0, An) are the associated 
group parameters. 

In order to describe the 2HDM potential, we introduce the 4- vector R'^: 



(2.15) 



a 



where S'^ in the full 8-dimensional field space must have the form: S'^ 
as required by SU(2)l gauge invariance. Moreover, as shown explicitly in Appendix [B], the 
imposition of U(1)y invariance and a Majorana constraint to be discussed below further 
reduces the form of the 4-vector matrices S'^ to 




(2.16) 



where O2 is the two-by-two null matrix. Consequently, in the Majorana scalar-field formalism. 
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we obtain for U(1)y- invariant 2HDM potentials that 

R'^ = R'^ . (2.17) 

However, we should stress here that if the U(1)y symmetry is lifted from the 2HDM potential, 
the 4-vector R'^ needs to be promoted to a 6-vector R"^ (with A = 0, 1, . . . , 5) and the 
corresponding structure of becomes non-trivial. In this respect, the Majorana scalar-field 
formalism has the analytical advantage in expressing the scalar potential of an U(l)- violating 
2HDM via a similar quadratic form with respect to R"^ as in (12. 5p for R'^. We discuss and 
demonstrate this application further in Section [61 

Under charge conjugation, the multiplet $ exhibits the following property: 

$ = C$* , (2.18) 

where C = ® cr° cr^, with C = C~^. Hence, $ satisfies a Majorana constraint, very 
analogous to the one obeyed by Majorana fermions. For this reason, we call this formalism 
the Majorana scalar-field formalism. In addition, the Majorana multiplet $ transforms under 
the reparameterization group GL(8,C) as 

$' = M$ , (2.19) 

with M G GL(8,C). However, as we will see below, the form of M cannot be general, but 
it is constrained by three basic conditions: (i) the conservation of SU(2)l symmetry by the 
transformation matrices M; (ii) the Majorana condition (I2.18P for any GL(8, C)-transformed 
multiplet (iii) the conservation of U(1)y symmetry by the transformation matrices M. 
Applying these three constraints on M, the 4-vector matrix S'^ is found to transform as 

e'/^A^S^ = M'^S'^M, (2.20) 

implying that R'^ transforms into 

R'/^ = e'/^A'^R'' , (2.21) 

where = det [M+M] > and A't G S0(1, 3). 

Since M G GL(8,C), the matrix M can then be represented in the full 8-dimensional 
scalar- field basis $ by the triple tensor product: 

M = M^^x a^'^a" . (2.22) 

As was mentioned in the introduction, there are two types of GL(8, C) transformations M 
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acting on $. The first one is a HF transformation, wliere tlie transformed multiplet 
transforms in the same way under SU(2)l as $, whereas the second one is a CP transformation 
where $' transforms in the same way as the charge-conjugated multiplet $*. Thus, for a HF 
transformation compatible with SU(2)l gauge invariance, we must have that M = U[ M Ul, 
where Ul is given in (12 .14^ . Instead, for a general CP and SU(2)L-invariant transformation, 
we must demand that M = Ul MUl- Consequently, the SU(2)L-invariant tensorial forms for 
the two types of transformation, which we denote as M±, are 



HF : M+ 

CP : M_ 



(2.23a) 
(2.23b) 



where we have used that V""" zcr^ V = za^, for any V G SU(2). 

It is now interesting to discuss the remaining two constraints imposed on the above 
SU(2)L-invariant structure of M±, resulting from the Majorana condition fl2.18p and the con- 
servation of the U(1)y hypercharge symmetry. The requirement that the Majorana condi- 
tion (12.181) should consistently hold for the multiplet $ and the HF/CP-transformed multiplet 
$' = M-i- $ produces the non-trivial constraint: 



= C M± C . 



(2.24) 



This last constraint reduces the form of the tensor M^j, defined in (I2.23ap and (]2.23b|) to 



M, 



( Moo Moi iMo2 Mo3 \ 

zMio iMu Mi2 zMi3 

iM2o iM2i M22 iM23 

\ iM^o iM^i M32 2M33 J 



(2.25) 



where all the components Mqo, Mqi, M02, . . . , M33 are real numbers. More details of this 
calculation are given in Appendix IB. 31 Thus, we observe that the Majorana condition applied 
to M reduces the reparameterization group from GL(8,C) to two subgroups isomorphic to 
GL(4,M), acting on a complex vector space. 

The HF and CP transformation matrices M-t should also respect the U(1)y hypercharge 
symmetry of the theory. Following a similar line of steps as for the SU(2)L-gauge invariance 
case, we require that M+ = UyM+UY, for a HF transformation, and M_ = UyM_ Uy, for 
a general CP transformation, where 



exp 



z6'y/2 (cr^ (g) a 



a 



exp {lOYa^/i^ ®(J^ ®(J^ , (2.26) 
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in the GL(8,C) representation, with 9y G [0,47r). Evidently, the above two constraints from 
requiring U(1)y invariance result in the commutator and anti-commutator conditions 



{M_ 



a 





(2.27a) 
(2.27b) 



for the HF and CP transformations, respectively. Since commu- 
tator relation fl2.27ap becomes M^^ [a'^, a^]® a" ® = It is not difficult to see that only 



= 0, 3 satisfy the last commutator relation, whereas Mic 
Then, M^y takes on the form: 



= 0, for a 



0,1,2,3. 



M, 



Mo3 




^M33 / 

leading to the following structure for the HF transformation matrix M+: 



00 




^M3o 



Moi 





^Mo2 




M32 



(2.28) 




(2.29) 



where 



Moo + Mo3 + iMso + iM33 Moi + M02 + iM^i - iMs2 
Moi - M02 + zM3i + ZM32 Moo - Mo3 + ^M3o - ZM33 



(2.30) 



is a general complex 2x2 matrix. The matrix form fl2.29p for M+ also provides closure in 
the 4- vector space of R'*, through the relation: 



(A- 



(2.31) 



where 6^^+ = det [T;T+] > and (A+)''^ G S0(l,3). 

Correspondingly, the anti-commutator relation given in ( ]2.27bl) leads to the constraint: 
Mfj,i, {a^, a^} ® a'^ ® (~^cr^) = 0. One can readily observe that only yU = 1, 2 satisfy the last 
anti-commutation relation, whilst Moa = M3Q = 0, for a = 0, 1, 2, 3. Thus, M^;^ acquires the 
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First Transformation Type 


Second Transformation Type 


Composite Type 


M+ 


M+ 


M+ 


M+ 


M_ 


M_ 


M_ 


M+ 


M_ 


M_ 


M_ 


M+ 



Table 1: Transformation properties after two successive operations ofM±. 



form: 

















\ 




iMio 


iMu 


Mi2 


zMi3 






iM2o 




M22 


IM23 




v 














/ 



(2.32) 



The resulting matrix M_ for general CP transformations is given by 




(2.33) 



where T_ is a complex two-by-two matrix given by 

/ Mio + Mi3 - iM2o - iM23 Mil - Mu - zMsi + iM22 \ 
T_ = . (2.34) 

\^ Mil + M12 - 1M21 - 1M22 Mio - Mi3 - zMso + iM23 J 

As before, the block-off diagonal form of M_ provides closure in the 4-vector space of R'^, 
since 

mLe^M_ = e"-/®(A_)tS^ , (2.35) 

with e"- = det [T1T_] > and (A_)''^ G S0(l,3). 

In addition we note that mixed transformations involving both M_|_ and M_ do not 
provide closure within the 4-vector space of R'^, i.e. 

M:^S'^M_ 9^ A'^S'^ . (2.36) 

Hence, two distinct S0(l,3) spaces exist which are compatible with U(1)y invariance. We 
denote these by (A+)'^ and (A_)'^, and their respective field transformation matrices by M+ 
and M_. Of course, combined transformation of different types are also possible, resulting 
in a composite transformation described by M+ or M_, as shown in Tabled! 

In summary, the HF and CP transformation matrices M± may be written down in the 
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HF/CP Symmetry 



Transformation Matrix T+ 
in the Basis (0i , ^2) 



Transformation Matrix T_ 
in the Basis (0i , 02) 



1 
1 



1 
-1 



U(l) 



PQ 



e-*" 


a e [0,7r) 



SU(2)hf/z, -S0(3) 



HF 



e cos 6* e sin 6 

— e*'^ sin 6 e*" cos 6 
e,a,P e [0,7r) 



CPl 



1 
1 



1 
1 



CP2 



1 
1 



1 
-1 



CPS 



COS 9 sin 9 

— sin 9 cos 9 

9 e [0,7r) 



cos 9 sin 9 

— sin 9 cos 9 
9 e [0,7r) 



Table 2: Matrix representations ofT± for 6 generic HF and CP symmetries of the 2HDM. 



following tensorial forms 
HF: M+ = 
CP : M_ = 



(a° + a') 
2 

(a' + la^ 



»T+ + 
®T_ - 



t; 



2 2 
Given the above representation of the HF and CP transformations, we observe that 



(2.37a) 
(2.37b) 



M^ 



CM 



(2.38) 



provided we set T_ = T^. This means that a general CP transformation can be thought of 
as a combination of a HF and a standard CP transformation. This is also consistent with 
the geometric interpretation presented in [25]. Likewise, the action of two successive CP 
transformations is equivalent to a single HF transformation, as can be seen from the last line 
of Table H 

In Table [21 we display the matrix representations of T+ (T-t) for the HF (CP) symme- 
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tries that we will be analyzing. In detail, the HF transformation matrices T+ are displayed 
in the second column of Table |21 These are the Z2 discrete symmetry [31], the Peccei-Quinn 
symmetry U(l)pQ [3S] and the HF symmetry S0(3)jjp [3T1 [221 IMl 125] which is isomorphic 
to a SU(2)jjp/Z2 transformation of 0i^2- Table |2] also exhibits the transformation matri- 
ces T± for three typical CP symmetries of the 2HDM potential: the CPl symmetry which is 
equivalent to the standard CP transformation 0i(2) — )■ 0*^2) [201 Ell EE], the CP2 symmetry 
where 0i(2) — )■ (~")02(i) ^"^^ the CPS symmetry which is a combination of CPl with an 
SO(2)hp/Z2 transformation of the Higgs doublets (j)i,2 [221I231IMI25] . 

Let us comment on the domains of the group parameters shown in Table |21 Specifically, 
we have considered a G [0, vr) for the U(1)pq symmetry, 6 G [0, tt) for the CPS symmetry, 
and a, (3,9 G [0,7r) for the SO(S)jjp symmetry. The parameter intervals for the potential 
symmetry groups are chosen so as to avoid double covers of the total symmetry group G, 
because of the presence of the SM gauge group SU(2)l ® U(1)y, and especially of U(1)y 
hyper charge [25] . 

Another important comment is in order here; for each CP symmetry, there should be a 
HF symmetry associated to it. This arises when the CP symmetry is raised to even powers and 
guarantees closure of the symmetry group (cf. Tabled]). For the CPl and CP2 symmetries, an 
even number of applications of the symmetry results in the identity mapping, i.e. (CPl)^"^ = I 
and (CP2)^'^ = I. However, for CPS, we obtain a non-trivial HF symmetry, i.e. (CPS)^"^ = 
SO(2)jjp/Z2. Unlike the CP symmetries, HF symmetries close within themselves, as shown 
in Table [H In Section 12. SI we will discuss further theoretical issues related to the breaking of 
the symmetry group G into a subgroup H, as these issues are important in order to generate 
the entire vacuum manifold associated to a given 2HDM potential. 

If the 2HDM potential V is invariant under a particular HF or CP symmetry Ghf/cpj 
realized by the matrices {A±)^^, then the theoretical parameters and L^;^ satisfy the 
relations: 

= M^At, (2.S9a) 
L,^ = L^uA^^A"^ , (2.S9b) 

Here, for convenience, we drop the subscript ± from (A-t)'^ and have implicitly assumed that 
0" = or e'^/^ = 1. Hence, for each HF or CP transformation acting on the Majorana field 
multiplet there is an equivalent transformation on R'^, as given in fl2.2ip . The tensor A'^ 
in the S0(1, S) space has then the following matrix form: 

Ahf/cp = diag (^1 , Ohf/cp) , (2.40) 
where Ohf/cp is a subgroup of 0(S) for the HF and CP symmetries under consideration. 
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HF/CP Symmetry 



S0(3) 



HF 



CPl 



CP2 



CP3 



rj 
^2 




/ 1 \ 

1 




/ -1 \ 

0-10 




U(1)PQ 




I 1 J 

a G [0,7r) 





Chf/cp Matrices in the Basis (R^ , , 



C2aCg — C2/3Sg 



'S2aCe — S2j3sl 



■S26»Ca-/3 —S2eSa-(3 

e,a,/3 e [0,7r) 



'S2eCa+l3 
'S2eSa+l3 
029 



1 
1 
1 



1 
0-10 
1 



1 \ 
1 



^ C28 

1 

\ -S20 C2e J 
a 



( -1 0\ 

0-10 
\^ 0-1 



. v - 

Q G [0,7r) 



/ C20 S2e \ 
0-10 
S2e C20 y 



Table 3: Matrix representation o/(9hf/cp defined in (12.401) for the 6 generic HF/CP symme- 
tries of the 2HDM potential. Here we use the shorthand notation cq = cos 6 and sq = sm6. 



In Table [3], we give the matrix representation of Ohf/cp; for the three HF and the three CP 
symmetries, respectively. 
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2.3 The Vacuum Manifold 



After minimization of the 2HDM potential, the field multiplet $ acquires, in general, a 



non-zero VEV, i.e. 



( 01 \ 



2^* 



la 



V2 



(2.41) 



where Vi^2 denote the VEVs of the Higgs doublets 0i,2- Employing the freedom of the 
SU(2)l ® U(1)y gauge transformations, the VEVs Vi,2 can be parameterized as: 



V2 = ^ 




(2.42a) 
(2.42b) 



where the vacuum manifold parameters fg, and ^ are all real. This parameterization of 
Vi 2 represents a single point of the vacuum manifold in the $-space, which we denote as fo- 
under this particular parameterization of the VEVs of the two doublets 0i,2, the equivalent 
extremal point in the R'^ basis in terms of the vacuum manifold parameters is: 



V1V2 cos ^ 
f °f 2 sin ^ 



(2.43) 



Our aim is to determine the entire vacuum manifold tM^ of the 2HDM potential, which 
amounts to finding all topologically distinct points of $, by appropriately acting on $0 with 
the set Ai that leaves the minimum of the 2HDM potential Vmin invariant. Thus, our task 
is to find Ai and its topological properties. We are interested in neutral vacuum solutions 
where both VEVs of 0i,2 are non-zero, i.e. situations where the vacuum component in 
fl2.42bp vanishes, = 0, and both 7^ 0- As a consequence of the latter, the VEVs Vi_2 
are invariant under rotations generated by the electromagnetic operator Qem = + Z/^ cr°, 
since QemVi^2 = (0,0)""", where = 1/2 is the hypercharge of 0i^2- Hence, if no HP or CP 
symmetries are present in the 2HDM potential, a non-trivial transformation of the VEVs 
Vi 2 can only be obtained by the action of the coset set: SU(2)l (8 \J{1)y/\J{1)^^. 
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If there is a HF (CP) symmetry group Ghf (Gcp) acting on the scalar potential V, then 
one needs to know whether there is a residual HF (CP) symmetry, Hhf (Hcp) say, which 
survives after spontaneous symmetry breaking. In such a breaking pattern: Ghf/cp ~^ 
Hhf/cp, the vacuum manifold point $0 is invariant under the action of the little group 
Hhf/cP) such that 

Hhf/cp : % ^ % = = % , (2.44) 

or equivalently Rg is invariant under Hhf/cp ? i-e. 

Hhf/cp: K -> K = (A")^.Ro = , (2-45) 

where Mh [(A^)^] is a representation of the unbroken group Hhf/cp in the GL(8, C) [S0(l,3)] 
space. As we will see in the next section, this is the case for the S0(3)hf model which breaks 
into the subgroup S0(2)hf — U(l)pQ. 

Consequently, a non-trivial HF/CP transformation of $0 or Rq can only be performed in 
the coset spaces: Ghf/Hhf or Gcp/Hcp. In a group-theoretic language, the vacuum manifold 
points $0 or Rq satisfying f l2.44p and f l2.45p are called orbit stablizers and the entire vacuum 
manifold can be generated by the transitive action of the total symmetry group G on them, 
where G = SU(2)l ® U(1)y ® Ghf/cp B Thus, in the GL(8, C) space, the entire vacuum 
manifold for a potential with HF / CP symmetry may be described by the set 

_^HF/CP ^ 1^ . $^_y^$^^ 7W G (GhF/Cp/HhF/Cp)® (SU(2)L®U(1)Y/U(l)em)} , 

(2.46) 

where $0 is the orbit stabilizer which is invariant under the little group U(l)cm ® Hhf/cp- 
The topological properties of A^^ or its generating set Ai under its homotopy groups, 
H„(A^), determines the nature of the topological defects [HI |12]. In particular, we have 
the existence of domain walls for Ho(A^) 7^ I, string solutions for Hi(A^) 7^ I, monopoles if 
Il2{M) 7^ I and textures if H„>2(A^) 7^ I [HI Il3], where I is the identity element. 

It is therefore vital to determine the representation of Ai in the full 8-dimensional $- 
space, for a HF and a CP symmetry. With this aim, we first note that a general element U 
of the SU(2)l ® U(1)y gauge group can always be written down as 

U = UlUy = expl^iOy^^ ®(T°®exp (^i^ly +2^2y^ cxp^^z^^y^ . (2.47) 

where Ul and Uy are given in fl2.14p and fl2.26p . respectively. Here, we also used the so-called 
Baker-Campbell-Haussdorf formula to factor out the third rotation due to the generator 



"^Throughout our study, we ignore the SU(3)^ colour gauge group which remains unbroken by the colour 
singlet VEVs of the Higgs doublets (f>i.2- 



17 



of Ul, where the transformed group parameters 6^'^'^ take values in the domain [0,47r). 
Using fl2.47p . one one can show that an element U-*- of the coset space SU(2)l ® U(1)Y/U(l)em 
may be represented in the $-space as 



a' -a' 



(g) 0-° (g) U_ 



with 



U± = exp ( i 0^ y + z ^2 y ) exp 



±i 



(2.48) 



(2.49) 



Note that the elements U"*" represent gauge transformations of the VEVs Vi_2 orthogonal to 
the U(l)gj^ electromagnetic group. In the $-space, the latter group can be represented by an 
expression very analogous to fl2.48p . where the 2x2 matrices U± are replaced with 



exp 



±i 



ey + e'- 



(2.50) 



Obviously, U"*" does not account for redundant rotations within U(l)gjjj, since °" Vi^2 = 
(0 , O)""" and ia'^Y\ 2 = (0 5 ^Y' ■ ^^is decomposition of the electroweak gauge group 
SU(2)l ® U(1)y into the electromagnetic group 11(1)^^^^ and the coset space U"*", the linear 
combinations 6± = ^{9y ± 9^) should be regarded as independent parameters which assume 
values in the domain 9± G [0, 27i). 

Given the representation f l2.48p for U-*", a non-trivial HF and CP transformation of the 
vacuum manifold point $0 is given by the GL(8, C) matrices 



Mi 



r ® u+ + 



0"" 



r* ® u. 



(2.51a) 



T- ® 



-102) 



(2.51b) 



where 7+ G Ghf/Hrf and 7± G Gcp/Hcp, with 7± being 2x2 complex matrices. Similarly, 
G Ghf/Hhf and G Gcp/Hcp are GL(8,C) matrices acting on the HF/CP coset 
spaces, whose tensorial form is very analogous to those given in fl2.37ap and fl2.37bp . 

At this point, it is important to reiterate that a HF symmetry Ghf of the 2HDM 
potential is closed under HF transformations M+ only, whereas a CP symmetry requires 
both types of HF and CP transformations M-t in order to obtain group closure, according to 
Table [TJ Likewise, the entire vacuum manifold for a 2HDM potential with a HF symmetry 
can be generated by acting only with transformation matrices of type M.+ given in fl2.51ap on 
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the initial vacuum manifold point $o- Instead, for a general CP-symmetric 2HDM potential, 
the complete vacuum manifold requires the use of both types of transformation matrices Ai± 
acting on $0 [cf. fl^naj) and fl2.51bp ]. 

As was already mentioned above, we may obtain an alternative description of the vac- 
uum manifold in the space. In this bilinear field basis, the entire vacuum manifold can 
be generated by the transitive action of the full group G on a single vacuum manifold point 
Rq, which is invariant under the orbit stabilizer group Hhf/cp [cf- f !2.45p ]. For this purpose, 
we would need to use the Ahf/cp or Ohf/cp matrices presented in Table [3] associated with a 
given HF/CP symmetry of the 2HDM potential. The vacuum manifold is then given by the 
set 

MfJ'''' = {R' : R'^ = AtRo> At e Ahf/cp/AIf/cp} , (2-52) 

where A^p^^^p is a possible residual HF/CP symmetry that remains intact after spontaneous 
symmetry breaking. In the gauge-invariant bilinear field basis, the SU(2)l ^ U(1)y gauge- 
group rotations are not present, so the nature of the topological defect solution depends only 
on the homotopic group properties of the coset bilinear field spaces: Ahf/A^f or Acp/A^p. 
We have checked that the analysis of the homotopy groups of the vacuum manifolds in the 
Majorana- field and the bilinear- field bases, A^^^^^^ and Ai^^^^ , lead to identical results. 

Finally, we should note that the breaking of the SM gauge group SU(2)l U(1)y to 
U(^)cm Siv6s rise to a vacuum manifold, which is homeomorphic to S^. This would imply 
that 113 [S"^] = Z, which would be indicative for the formation of non-trivial topological 
configurations called textures. However, such local textures turn out to be gauge artifacts 
since they can be removed by a gauge transformation [H]. Global textures and monopoles, 
whilst unstable due to Derrick's theorem, can be cosmologically interesting, for instance 
global monopoles can provide a mechanism for structure formation [HU] . For this reason, our 
focus will be on non-trivial topological configurations that arise from the breaking of HF or 
CP symmetries: Ghf/cp Hhf/cp- 



3 Neutral Vacuum Solutions of the HF Symmetries 

We start our analysis by considering the three generic HF symmetries: Z2, U(1)pq and 
S0(3)jjp. These HF symmetries impose specific relations [23] among the parameters of the 
2HDM potential, which are presented in Table |H For the Z2 symmetry, the quartic coupling 
A5 can always be made real by a simple phase redefinition of one of the two Higgs doublets 

01,2- 

Given the constraints on the potential parameters due to the HF symmetries, the four 
general convexity conditions (I2.10ap - (l2.10cl) and (lC.6ip become greatly simplified. These four 
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Symmetry 






Ai A2 


A3 A4 A5 


Ae 


A7 


Z2 









- - Real 








U(1)PQ 









- - 








S0(3)hf 







- Ai 


- 2Ai - A3 









Table 4: Parameter relations in the 2HDM potential that result from the imposition of the 
three generic HF symmetries. A dash indicates the absence of a constraint. 



Convexity Condition 


Z2 


U(1)PQ 


S0(3)hf 


1 


Ai > 


Ai > 


2Ai > IA3I 


2 


A2 > 


A2 > 




3 


2VA1A2 > lAsI 


2VA1A2 > A3 




4 


A4 > A5 


A4 > 





Table 5: The four convexity conditions for a hounded- from-helow 2HDM potential for each of 
the three HF symmetries. A dash signifies the absence of any additional constraints on the 
parameters. 

conditions are exhibited in Table [51 In the S0(3)jjp case, the convexity conditions are not 
independent of each other and only one distinct condition survives. 

We will now derive analytical expressions for the neutral VEVs of 0i 2 for each of the 
three HF symmetries, by utilizing the Lagrange multiplier method. These results will enable 
us to study in more detail possible topological defects that can emerge from a non-trivial 
vacuum topology of the theory, as shown in Section [51 

3.1 Z2 Symmetry 

The discrete Z2 symmetry of the 2HDM is defined by the following transformations of the 
two Higgs doublets 0i_2: 

01 01 = 01 , 

02 02 = -02 • 

To solve the extremization condition fl2.9al) . we consider two cases: (i) det[N^j,] 7^ and (ii) 
det[N^,^] = 0. In the first case, the matrix N^j, can be inverted and the 4- vector R'^ can 
be straightforwardly derived, whereas in the second case N^j^ is not invertible and a slightly 
different strategy needs to be deployed to determine R^. 

Taking into account the parameter restrictions of Table [H for the Z2 symmetry, we may 
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now calculate the determinant of N^,^ (see also Appendix [C]) . This can be expressed in the 
factorized form: 

det[N^.] = [A2-(A4 + C)'] [(A3-C)'-4AiA2] . (3.1) 

For the Z2 symmetry, the extremization condition N^j.R'^ = decomposes into two separate 
matrix equations: 



Ai + A2 + A3 - C Ai - A2 \ / R° \ I 1^1 + 1^1 

Ai - A2 Ai + A2 - A3 + C / V / V Ai? - At2 

A4 + A5 + C 

A4 - As + C 




(3.2a) 



(3.2b) 



Assuming that N^;^ is non-singular, the above matrix relations can be inverted and the 
individual components of R'^ for an arbitrary point on the vacuum manifold are found to be 

= 4A,A. - (A3 - C)^ ' ^'-'"^ 

R^ = , (3.3b) 

R^ = 0, (3.3c) 

p3 2A2/x?-2Ai/ii + (A3-C)(/i?-/ii) „ ,^ 

= 4A,A2 - (A3 - 0^ • ^'-'^^ 

From the defining equation (12. 3p for the 4- vector R^, the following analytical expressions for 
the VEVs of the Higgs field bilinears are easily obtained: 

^'^^'^^^ = 4A,A2 - (A3 - C)^ - ° ' ^'-'"^ 

/A^A\ 2Ai/i2 - (A3 - C)Ai? rn A-u\ 

^'^^'^^^ = 4A,A2-(A3-C)^ - ° ' ^'-'^^ 

(0l</)2) = {4<Pi) = . (3.4c) 



In order to have a neutral vacuum solution, we must satisfy the condition f l2.12p . namely 
that R'* is a null 4-vector, with R'^R^ = 0. This restriction leads to 

4A1A2 - (A3 - 0^ " ' ^ ^ ^ 

which completely specifies the Lagrange multiplier. More explicitly, requiring that the nu- 
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VEV parameter 


Ci 


C2 


n 












V A2 



















Table 6: The two neutral vacuum solutions to the Z2 symmetric 2HDM potential for 
det[N^!y] 7^ 0. The Lagrange multipliers Ci,2 o'^e given in f l3.6ap and (]3.6b|) . 



merator of (13.51) vanishes, we find two solutions for the Lagrange multiplier: 

Ci = A3 - ^ , (3.6a) 
C2 = A3 - ^ . (3.6b) 

A*2 



Using the specific parameterization ( I2.42ap and (12.42bp for the VEVs of 0i^2) we can deter 



mine the vacuum manifold parameters (f°, f2, tij, ^) for the two values ^1,2 of the Lagrange 
multiplier given in f l3.6ap and f l3.6bp . The results are given in Table [61 Moreover, we have 
verified that the two solutions Ci,2 do not lead to a singular matrix N^,y. 

In order for a set of neutral vacuum solutions to correspond to a local minimum of the 
potential, we require that the Hessian of the Z2 invariant 2HDM potential is positive definite. 
The general Hessian of the Z2 invariant 2HDM potential with respect to v\ and is given 
by 

/ + 3^^(^0)2 + 1A345(^;0)2 A345^^M \ 

H = . (3.7) 

\^ A345« -l4 + ^Hvir + |A345(^^?)^ ; 

Here we introduce the common summation conventions between the quartic couplings of the 
model: Xab = Aa + A;, and Xabc = Aq + Af, + Ac. Thus, the positivity of H leads to additional 
constraints, which are listed in Table [71 Specifically, the first condition in Table [7| corresponds 
to having a local minimum, whilst the second one is to ensure that this minimum is the lowest 
one. If fif = ^2 Ai = A2, the global minimum is given by 

Vo = - ^ . (3.8) 
As can be seen from Table [6l when the determinant of N^;^ is non-zero, at least one of 
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Condition 


Ci 


C2 


1 


/i^ > 




2 


fii ^ 2Ai 

-^345 


A'l ^ -^345 
/i^ 2A2 



Table 7: Minimization conditions for two neutral vacuum solutions in a Z2 symmetric 2HDM 
potential, with det[N^i,] 7^ 0. The first condition corresponds to having a local minimum and 
the second one is for this minimum to be the lowest. 



the VEVs of the Higgs doublets 0i^2 must be zero, in order to have a neutral vacuum solution. 
As we will discuss in Section I3.1.2[ such solutions do not lead to topological defects, such as 
domain walls in this case, and they are not of interest for the present study. We now turn 
our attention to the neutral vacuum solutions that can occur when the matrix N^i, becomes 
singular for a specific choice of the Lagrange multiplier. 



3.1.1 Neutral Vacuum Solutions from a Singular Matrix N 

We now consider the possibility that the matrix N^,y has no inverse, by requiring that its 
determinant given in (13. ip vanishes. Equating separately the two factors in (13.11) to zero, we 
obtain four solutions: 

Ci,± = -A4±A5, (3.9a) 



C2,± = ±2^/xjr2 + X^. (3.9b) 

Since the extremization condition for the Z2 invariant potential splits into two separate 
matrix equations, (I3.2ap and f l3.2bp . the application of either of the above four Lagrange 



multipliers only results in one of the matrices in the equations becoming singular. For the 
solution C2,±; it is the 2x2 matrix in (I3.2ap which becomes singular. However, since the 
RHS of (I3.2ap is in general a non-zero vector in this case, unless /^i = /^i = 0; this matrix 
equation is overdetermined. Unless the parameters /if 2 the quartic couplings Ai,2,3 satisfy 
an unnatural fine-tuning relation, the matrix equation (I3.2ap becomes incompatible for the 
Lagrange multiplier C2,±- We therefore reject the second solution (^2, ± and focus on the first 
solution Ci,±- 

For the Lagrange multiplier solution Ci,±; the matrix in f l3.2bp becomes singular, whilst 



the matrix equation (I3.2ap can be inverted in general, using standard linear algebra methods. 
Evaluating the singular matrix in f l3.2bp . we observe that the solution yields = 0, but 



leaves undetermined. Likewise, the solution (^1 _ renders R^ = 0, but R^ 7^ in general. 
The two solutions are related by a reparameterization of the doublets, since $2 implies 
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— )• Ci,-- Therefore, only one solution of the Lagrange multipliers needs to be considered. 

Having the above in mind, we consider the solution Ci,-; where A5 enters additively in 
all resulting equations. Substituting ^1 _ into f l3.2ap gives 



4A1A2 - A345 

3 2A2Aif - 2Ai/i2 + A345(/if - fil) 

K = 72 • {6.10b) 



4A1A2 - A345 

In terms of field bilinear VEVs, R° and imply that 



t , \ 2A2/ii - A345/i2 „ fQTT \ 

= , To — > , (3.11a) 



4A1A2 - A 



2 
345 



mh) = —TT-^ 72 — > 0. (3.11b) 

4A1A2 - A345 

In addition, the constraint R^ = translates into {(p\(p2) = {4>\'Pi)^ which can only be satisfied 
if the phase ^ is a multiple of vr, i.e. ^ = nvr, with n being an integer. 

In order to uniquely fix the undetermined component R^, we require now that R^ is a 
null vector, i.e. R^R'^ = 0. Employing this last condition, we find that 



/pl\2 ^ 4 [2A2/ii — A345/i2] [2Ai/i2 — A345/i^] 

[4A1A2 - Ai 1' 



'345J 



Comparing fl3.10ap . (]3.10bp and fl3.12p with the R'^ parameterization in f l2.43p with 11^ = 
and = 0, we obtain 



4A2/i? - 


2A345/i2 


4A1A2 


^^345 


4Ai/ii - 


2A345/i^ 


4A1A2 


^345 



(3.13a) 
(3.13b) 

By analogy, we may calculate the vacuum manifold parameters related to the Lagrange 
multiplier Ci + = — A4 + A5. These are found simply by replacing A345 in all equations with 
A345, where we extended the summation convention as: Xabc = + \ ~ K- As we discuss in 
Section [531 the space of the entire vacuum manifold is generated via the transitive action of 
the total symmetry group on this particular set of the vacuum manifold parameters. We have 
also checked that the VEVs of the Higgs doublets 0i,2 obtained by the Lagrange multiplier 
method satisfy the extremization conditions given by the usual tadpole equations f l2.2ap 
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-1.0 -0.5 0.0 0.5 1.0 




Figure 1: Contour plot depicting the shape of the Z2 invariant 2HDM potential V for the 
parameter set {/if , yU^) -^i) -^2, A345} = {1,1,1,1,1}, in arbitrary mass units and normalized 
such that Vmin = 0. The four degenerate and disconnected global minima are shown in black 
around the central local maximum. The four minima form two pairs; the members within 
each pair are related by the Z2 symmetry and the two pairs are related to one another by 
U(1)y. 



and (IlJb]). 



To determine whether the above extremal solutions represent local minima as well, we 
require that the Hessian H in ( 13 .Tp . evaluated at the extremal points, is positive definite. 
This requirement generates two conditions: 



4A1A2 - A 



> . (3.14b) 



345 



These two inequalities are equivalent to the positivity conditions for the squared VEVs in 
(13.1 lap and (I3.11bp . provided 4A1A2 > AI45 and Ai > 0. The constraint Ai > represents one 
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of the convexity conditions for the Z2-symmetric 2HDM potential (see Table ED- However, the 
restriction 4A1A2 > A345 has not been accounted before and creates two additional inequalities 
from the numerators of the fractions given in (13. Hal) and fl3.11bl) . These can be summarized 
in the double inequality 

ZA2 /i2 ^345 

Comparing this double inequality with the second line in Table [TJ we see that local minima 
with f°2 = and f°2 7^ cannot coexist. The value of the potential at the local minimum 
associated with the Lagrange multiplier „ is given by 

= JTx ■ ^^-^^^ 

4A1A2 - A345 

The corresponding value Vq for the local minimum related to = —-^4 + -^5 is obtained 
by making the substitution A345 — t- A345 in (I3.16p . Between these two solutions, the lowest 
minimum is given by = — A4 + A5, if A5 > 0, and by (i- = — A4 — A5, if A5 < 0. Hence, 
the potential at the lowest minimum is given by 

y ^ (A3 + A4 - IA5I) nlfil - Ai/x^ - X2fit ,g 
° 4AiA2-(A3 + A4-|A5|)' 

Note that this lowest minimum becomes a global one of the Z2-symmetric 2HDM potential, 
if (I3.15P is fulfilled. Otherwise, the global minimum is given by (13. 8p . A numerical example 
of a Z2-symmetric 2HDM potential, where both are non-zero, is shown in Figure [H 



3.1.2 Z2 Topology 

It is now important to determine the topology of the vacuum manifold for the Z2 invari- 
ant 2HDM potential, applying some of the general results presented in Section 12.31 In the 
symmetric phase, the Z2 invariant 2HDM potential is governed by the total symmetry group 
Gza = Z2 ® SU(2)l U(1)y, including the electroweak gauge group. After spontaneous 
symmetry breaking of the electroweak gauge group, we have 

SU(2)l ® U(1)y ^ X S" ^ U(l),^ ~ S' . (3.18) 



In the above, we used the well-known homeomorphisms between compact groups and n- 
spheres denoted as S"" (or S""): U(l) ~ and SU(2) ~ S^. According to our discussion in 
Section [273| in the absence of any HF/CP symmetry in the theory, the vacuum manifold of 
the 2HDM will then be homeomorphic to the coset space {S'^ x S'^)/S^, which in turn is 
homeomorphic to S^, i.e. {S'^ x S'^)/S^ ~ S^. 
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In the present case, there exists an additional discrete Z2 symmetry acting on the 
2HDM, which can break to the identity, i.e. Z2 — )■ I, after electroweak symmetry breaking. 
If this happens, the breaking pattern of the total symmetry group proceeds as follows: 

Gz, = Z2®SU(2)l®U(1)y ^ Hz, = I®U(lU. (3.19) 

As a consequence, the topology of the vacuum manifold will then be described by the coset 
space Ml"" = Gza/Hzj. 

In order to generate the complete set A^|^ of the vacuum manifold points in the 
space, we first need to find an initial point ^0 of the Majorana scalar-field multiplet, which 
remains invariant under the little group Hz,. Then, A^|^ will be generated by the transitive 
action of Gz, on $o- In the parameterization of the Higgs-doublet VEVs Vi,2 of fl2.42ap and 
fl2.42b|) . the Majorana scalar-field vacuum point $0, which is invariant under Hz, = U(l)em' 
is given by f ^ = and ^ = 0. 

Let us first consider the non-trivial case where 2 7^ as discussed in Section 13.1.11 
The general vacuum manifold point $ is given by 

$ = Ml^ $0 , (3.20) 

where the HF transformation matrix Ai'^ is stated in fl2.51ap and 7+ = T_(_ = {o"" , a^} are 
the 2x2 HF transformation matrices given in Table [2] under the Z2 symmetry. It is interesting 
to see the different roles of the Z2 symmetry and the U(1)y hypercharge symmetry, according 
to the more intuitive chart: 



(::) 




(:::) 








Z2t 




|Z2 








V , 




. V2 ) 



Observe that for Z2-symmetric 2HDM scenarios with two non-zero VEVs f°2 7^ 0, we cannot 
move via a U(1)y transformation from one vacuum configuration, e.g. {v^ , f^), to its Z2- 
symmetric one, i.e. (f ° , — ^2) or (— f° , f^)- However, if or V2 were zero, then such a 
transformation would be possible, and the discrete vacua will be connected via a continuous 
U(1)y gauge transformation. In the latter case, there are no topological defects, such as 
domain walls or superconducting condensates similar to the ones discussed by Hodges jM] . 
even though such scenarios might be interesting as they predict stable scalars which may act 
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as DM (see, e.g. 00]). 

On the other hand, the Z2 invariant 2HDM, where the two VEVs are non-zero, can 
lead to non-trivial topological solutions, such as domain walls @. The vacuum manifold in 
the •I'-space may be given by 

Ml' ~ Z2xS\ (3.22) 

where the second factor comes from the breaking pattern of the electroweak gauge group 
as given in fl3.18p . Thus, the action of the zeroth homotopy group on this vacuum manifold 
is non-trivial, since Hq [Z2 x S^] = Uq [Z2] x Ho [S^] ^ I, with Hq [5^] = I p]. This leaves 
the possibility for the formation of domain walls in the Z2 symmetric 2HDM, whose spatial 
profile is studied in Section |5l 



3.2 U(l)pQ Symmetry 

We now analyze the Peccei-Quinn symmetry of the 2HDM, which is defined by the following 
transformations of the two Higgs doublets (f)i^2- 

^1 -> = e-^"(/.i , 
02 ^ 02 = e^"02, 

where a G [0, vr). The study of the neutral vacuum solutions of the U(1)pq invariant 2HDM 
proceeds in a very analogous fashion to the Z2 invariant 2HDM discussed in the previous 
section, since the only additional parameter restriction in the U(1)pq invariant theory is that 
one now has A5 = 0. Therefore, we only quote a few key results here. 

For neutral vacuum solutions resulting from a non-singular matrix N^,^, the VEVs are 
given by fl3.4ap and fl3.4cl) . with A5 = 0, i.e. 

= ^T— T 7T > , (3.23a) 

4A1A2 - (A3 - C)"' 

(0202 = -TT—. 7T > , (3.23b) 

4A1A2 - (A3 - C)"' 

(0I02) = (0^01) = . (3.23c) 

There are two Lagrange multiplier solutions for this situation which are given by fl3.6ap and 
(13.6bp . Because of this close similarity, the vacuum manifold parameters are exactly the same 



as those detailed in Table |6] of Section [3?T1 Correspondingly, the conditions for each solution 
to correspond to a minima are given in Table [3 As in the Z2 case, the U(l)pQ-invariant 



^Here we assume that there are no other sources that violate the Z2 symmetry of the theory, e.g., either 
by Yukawa coupUngs, or by anomahes |52j . 
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2HDM must also have at least one doublet with a zero VEV, when det[N^,y] 7^ 0, which only 
leads to topologically trivial configurations. We are, therefore, only interested in neutral 
vacuum solutions for which the matrix N^,^ is singular. 



3.2.1 Neutral Vacuum Solutions from a Singular Matrix N 

In order for the matrix N^i, to have no inverse in the case of the Peccei-Quinn symmetry, we 
require that the expression given in (13. ip with A5 = be equal to zero. This requirement 
leads to two candidate solutions: 

Ci = -A4, (3.24a) 



C2,± = ±2^/xJ^2 + Xs. (3.24b) 

However, for the same reasons as in the Z2 case, we have to reject the second solution (2,±, as 
it leads to an incompatible matrix equation, unless there is a particular fine-tuned relation 
between the parameters of the 2HDM. Therefore, we only focus on the first solution (i. 

Under this choice for the Lagrange multiplier both and remain undetermined, 
since the 2x2 matrix in ( I3.2bl) becomes the null matrix. The remaining components of the 



vector R^ are found using fl3.2al) and have the form: 

2Ai/i^ + 2Ai/i?-A34(/U? + /ii) 

— -T — ; 7-^ 



4A1A2 — A34 

2A2/i? - 2A2/i2 + A34(/U? - /U2) 



(3.25a) 



From these expressions, we obtain by means of (12.31) the VEVs of the scalar-field bilinears 

('^l'^^) = ^ ' (3.26a) 

4A1A2 — A34 

{<Pl4..) = ''f~'lf > . (3.26b) 

4A1A2 — A34 



For a neutral vacuum solution we require, as before, that R'^ satisfies R^R'^ = 0, which leads 
to the relation: 

/oU2 I m2N2 ^ 4 {2X2f4 - ^34/^1) (2A1/X2 - A34/if ) 

(4A1A2 - Ai^)^ 

Employing the parameterization of R'^ in (I2.43P , we find that the vacuum manifold parameters 
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for the Lagrange multiplier d with f ^ = are: 





2A34/i2 


4A1A2 




'4Ai/xi - 


2A34/if 



^ e [0,27r) . (3.28c) 

Notice that the phase ^ remains undetermined, signifying the presence of a massless Goldstone 
boson, the so-called PQ axion 

The conditions for a global minimum are identical to those of the Z2 case with A5 = 0. 
Thus, we have a global minimum with f?2 7^ 0) provided 

^^'^ < 4 < ^ • (3-29) 



2A2 /^2 ^ 
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The value of the U(l)pQ-invariant 2HDM potential at the global minimum is given by 



•^34/^1/^2 ~ -^1/^2 ~ -^2/^1 

4A1A2 — AI4 

As before, we find that neutral vacua where f 5* 2 7^ and f ° = or f 2 = cannot co-exist 



Vo = ^^"7, '^"^,2 • (3-30) 

4A1A2 — A34 



3.2.2 U(l)pQ Topology 

Let us now discuss the topology of the vacuum manifold associated with the U(l)pQ-invariant 
2HDM potential. The total symmetry group of the potential in the symmetric phase is 
Gu(i)pQ = U(l)pQ ® SU(2)l ® U(1)y- After electroweak symmetry breaking [cf. fl3.18p ]. the 
U(l)pQ symmetry breaks into the identity I, so the unbroken group is Hu(i)pQ = I (g) U(l)em- 
As a consequence, the vacuum manifold in the $-space is given by the set 

mT^'-'' = Gu(i)pq/Hu(i)pq ^ S'xS\ (3.31) 

where we used the fact that U(1)pq is homeomorphic to S^. We now observe that the first 
homotopy group of this vacuum manifold is non-trivial, i.e. Hi [S^ x S^] = Hi [S^] x Hi [S^] = 
Hi [S^] = Z 7^ I, since Hi [S^] = I. This implies that the U(l)pQ-invariant 2HDM has a string 
or vortex solution, which we analyze in Section |5l 

It is interesting to discuss the construction of the vacuum manifold in the $-space. As 
stated in f l2.51ap . a general point of the vacuum manifold is given 
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where $o is defined in terms of tlie non-zero VEVs v^2 given in fl3.28ap and f l3.28bp and by 



setting ^ = 0. Moreover, in tlie 8-dimensional Majorana $-space, tlie HF transformation 
matrix A^^''^^^'^ takes on tfie form: 



M 



U(1)PQ 



(g) exp 



2ia 



a' -a' 



exp 



-2ia 



_(_^^L/ ,6 ,6- — a) 



(3.32) 



Here, we have exphcitly displayed the dependence of the gauge-group factors U± on their 
group parameters and made use of the fact that the HF transformation matrix 7+ = e"*""^^ 
for the PQ symmetry may be written asT+ = e-^°e^°(''°-'^'). We may now re-define the group 
parameter ^_ as ^_ = 0_ — a G [0,27r), and so having the group parameter 2a G [0,2-7?) to 
span now the complete space of the U(1)pq group. Note that this result is identical to the 
one that would be obtained in the space, as can be easily deduced from Table [3l 



3.3 SO(3)hf Symmetry 

An interesting HF symmetry emerges from the invariance of the 2HDM potential under a 
naive SU(2)hf transformation of the Higgs fields, i.e. 

01 ^ (p[ = e-^°cose0i + e"^^ sin ^02, 

02 ^ 0'2 = -e^^ sin ^01 + e*" cos ^02. 

To avoid a double cover of the SU(2)hf group because of the presence of U(1)y hypercharge 
rotations, we have to restrict the group parameters 6, a, P to lie in the interval [0, vr). Hence, 
the actual HF symmetry is the coset group SU(2)hf/Z2 [23], which is isomorphic to S0(3)hf 
in the field-bilinear space. For this reason, this symmetry was called the S0(3)hf sym- 
metry. 

The parameters of the 2HDM potential under the S0(3)hf symmetry are restricted, as 
shown in Table HI In fact, most of the results can be easily recovered from the Z2 case in 
Section \3A] by making the replacements: A2 \i, A4 — )■ 2Ai — A3, /i2 — )■ fif and putting 
A5 = 0. Therefore, we will only report key intermediate results in this section. As before, we 
first assume that the inverse of N^,^ exists, with the determinant of N^^, given by 

det[N^,] = (2Ai + A3 - C) (2Ai - A3 + C)' • (3.33) 
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Because of the more restrictive nature of the S0(3)hf symmetry, the matrix equation N^^^R*^ = 
sphts into four separate equations: 

(2Ai + A3-C)R° = (3.34a) 

(2Ai-A3 + C)R' = 0, (3.34b) 

(2Ai-A3 + C)R' = 0, (3.34c) 

(2Ai-A3 + C)R^ = 0. (3.34d) 

On the basis of the above assumption that N^,^ is invertible, the components of the 4-vector 
R'^ are easily found to be 

rO = , (3.35a) 

2Ai + A3 - C 

Ri = r2 = r3 = . (3.35b) 

On the other hand, tlie constraint for a neutral vacuum solution requires that R^ is a null 
vector, satisfying R^R^ = 0. Since all the "spatial" components R^'^'^ vanish, so should the 
"time" component, i.e. R° = 0. This last result tells us that the Higgs doublets should have 
vanishing VEVs, i.e. — 0; leaving the electroweak gauge group unbroken. This is an 
unrealistic scenario and can only be obtained in the limit /x^ — )■ 0, or — )■ ±00. We will 
therefore investigate neutral vacuum solutions that result from a singular matrix N^,^. 

3.3.1 Neutral Vacuum Solutions from a Singular Matrix N 

From f l3.33p . we readily see that the following choices of the Lagrange multiplier render N^i, 
singular: 

Ci = -2Ai + A3 , (3.36a) 
C2 = 2A1 + A3. (3.36b) 

However, from fl3.34al) . we notice that the solution (2 implies either R° — )■ 00, or fif — )• 0, 
both of which lead to unrealistic scenarios of electroweak symmetry breaking. Therefore, we 
concentrate on the Lagrange multiplier solution (i. 

Considering the Lagrange multipher solution ^1, we obtain from fl3.34ap that 



R° = ^ . (3.37) 



Instead, from fl3:34bl) -f l3:34dl) . we see that all the "spatial components" R ' ' remain unde- 
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termined. The only constraint that can be placed upon the three "spatial" components of 
R'^ is the requirement of a neutral vacuum solution, R^R'^ = 0, which implies that 



(Ri)2 + (R2)2 + (R3)2 = A_ . (3.38) 

4Ai 



In terms of the vacuum manifold parameters and ^, (13. 370 and ( 13.38P are translated into 

- sine, = ^ cos^^, (3.39) 



where ^ G [0, 27r) and 6 G [0, vr) remain undetermined. The latter signifies the presence 
of two Goldstone bosons. Specifically, the one associated with the phase ^ is a CP-odd 
scalar, whereas the one related to the polar angle 6' is a 'CP-even' boson. This result can be 
cross-checked independently from the explicit analytical expressions presented in [16] for the 
general Higgs-boson mass matrices. The global minimum of the S0(3)HF-symmetric 2HDM 
potential is given by 

V„ = - i| . (3.40) 

Such a global minimum is always guaranteed, as long as /if is positive and the bounded-from- 
below condition, 2Ai > IA3I given in Table [5] is satisfied. 

3.3.2 SO(3)hf Topology 

It is interesting to analyze the topology of the vacuum manifold arising from the spontaneous 
symmetry breaking of an S0(3)HF-iiivariant 2IIDM potential. In the symmetric phase of the 
theory, the S0(3)HF-invariant 2IIDM potential has the symmetry, which is described by the 
group 



Gso(3)hp = (SU(2)hf/Z2)®SU(2)l®U(1)y = S0(3)hf ® SU(2)l ® U(1)y • (3.41) 

Using the results of the previous section, we see that out of the three generators of the 
SU(2)hf/Z2 group, one linear combination of generators, ((J° + (J^)/2 related to a residual HF 
symmetry, which we call U(l)jjp, remains unbroken after the electroweak symmetry breaking, 
resulting in the little group 

Hso(3)hp = U(1)hf®SU(2)l®U(1)y = S0(2)hp SU(2)l ® U(1)y . (3.42) 
Then, the vacuum manifold Ai^^^^^^ may be described by the product of spaces: 

^so(3)h, ^ Gso(3)h,/Hso(3)h, ^ S'xS', (3.43) 
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where the first factor S*^ is obtained using the known homeomorphism SO(3)hf/SO(2)hf 
^ S*^ and the second one is due to the breaking of the electroweak group to U(l)^^. 

SOf3) 

We observe that the second homotopy group of is non-trivial. More explicitly, 

U2 [S^ X S^] = U2 [S^] X U2 [S^] = II2 [S^] ^ I, since U2 [S^] = I. Consequently, spontaneous 
symmetry breaking of the S0(3)jjp-symmetric 2HDM can give rise to global monopoles. 

As with the previous HF symmetries, we are able to construct the entire vacuum mani- 
fold by the transitive action of the total group Gso(3)jjp stated in fl3.4ip on the vacuum point 
$0) which remains invariant under the little group Hso{3)jjp given in (13.42 p . An appropriate 
representation of $0 consistent with the latter property is given by the VEVs 



= , V2 = -= , (3.44) 





where we set 6' = = in f l3.39p . The general point $ on the vacuum manifold is then 

SO(3) 

given by the action of the coset set of HF transformation matrices on i-e. 

$ = M^^^^™ $0, where 

Here, the 2x2 HF transformation matrices 7+ belong to the coset space of the S0(3)hf sym- 
metry in the adjoint representation, i.e. 7+ G (SU(2)jjp/Z2)/U(1)hf, and can be represented 

as 

e~*" cos 6 e~*^ sin 6^ \ / e*'^ \ / cos 6 sin 6 

r, = { ] = e-^ 

-e^^sin^ e^"cos^ / \ 1 / \ -e*("+'^) sin ^ e^^^+'^^cos^ 

(3.46) 

where we set the free U(1)hf phase x to be x = « — /3, in obtaining the second equation. As 
in the PQ symmetry case, the overall factor e"*'' can be absorbed into the definition of the 
gauge-group parameter 6-, i.e. by defining 9^ = 9^ — (5 E [0,27r). The HF transformation 
matrices A^^^''^^^'" can then be written down as 

Mf^'^-- = I ^4^] ® ^\J^{9\9',9.) 

n-O — n-'^A / COS 6' sin^^ , _ _ 

+ ( — ^ ® U)U_(^^^^^_) . (3.47) 

gi^ g g-i(a+/3) ^qs 9 ^ 

If we ignore the gauge rotations by setting U-t = a", the action of J\4^^^^^p on $o then 
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Symmetry 




^2 


m^2 Ai 


^2 A3 


A4 A5 


Ae 


A7 


CPl 






Real - 




- Real 


Real 


Real 


CP2 









Ai - 






-As 


CPS 









Ai - 


- 2X1 — A3 — A4 









Table 8: Parameter relations in the 2HDM potential that result from the imposition of the 
three generic CP symmetries. A dash indicates the absence of a constraint. 



Convexity 
Condition 


CPl 


CP2 


CPS 


1 


Ai + A2 + A3 > 


2Ai > -A3 


2Ai > IA3I 


2 


A4> A5+^,+^,+^^ 


A4 > — i?5 


2A4 > 2Ai — A3 


3 


A4 > A5 


A4 > A5 ^ 




4 


A1A2 — 4AI > 
A1A2+A2A2-A3A6A7 


2Ai - A3 > 

i\X6\^{XA-Ri)-&h{hR<i-Rs.h) 




A4+A5 


AMAsP 



Table 9: The four convexity conditions for a hounded- from-helow 2HDM potential for each of 
the three CP symmetries. A dash signifies the absence of any additional constraint. 

generates the general vacuum manifold point given in (13.391) . with ^ = a + /3 G [0,27r) and 
9 G [0, tt). Thus, the vacuum manifold of the SO(3)HF-broken 2HDM is homeomorphic to S"^, 
parameterized by the azimuthal angle 9 and the polar angle ^ = a + This parameterization 
will be used in Section |5] to analyze the monopole solution in this model. 



4 Neutral Vacuum Solutions of the CP Symmetries 

In this section, we will study the three generic CP symmetries, termed CPl, CP2 and CP3. 
These three CP symmetries impose specific relations [23] among the parameters of the 2HDM 
potential, which are presented in Table 

Implementing the constraints on the potential parameters due to the CP symmetries, 
the four general convexity conditions fl2.10ap - fl2.10cp and flC.6ip take on a simpler form. These 
four conditions are displayed in Table [91 In particular, for the CP3 case, the four convexity 
conditions are not all independent of each other, so only the two distinct conditions are 
presented. 

As in the previous section, our aim is to derive analytical expressions for the neutral 
VEVs of 01^2 in terms of the 2IIDM potential parameters for each of the three CP symmetries, 
by making use of the Lagrange multiplier method. These results will be used to determine 
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the existence and the nature of possible topological defects which will be discussed in detail 
in Section Ei 



4.1 CPl Symmetry 

The discrete CPl symmetry of the 2HDM represents the standard CP transformation of the 
two Higgs doublets 0i,2, given by 



t>2 -> 



'1 ) 



Taking into account the CPl parameter restrictions of Table [HI we calculate the VEVs of 0i^2 
by imposing the extremization condition (I2.9ap and the condition (I2.9b[l for an electrically 
neutral vacuum. As before, we consider two cases: (i) det[N^,^] 7^ and (ii) det[N^,^] = 0. 

The determinant of N^j, resulting from a CPl-invariant 2HDM potential follows from 
Appendix [C] and can be expressed in the factorized form: 

det[N^,] = (A45 + C)[(A45 + C)(4AiA2-(A3-C)')-4AiA2-4A2A? + 4A6A7(A3-C)] ■ 

(4.1) 

Moreover, the extremization condition N^j^R'^ = decomposes into two equations: 



^ A1 + A2 + A3 — Ag + Ay Ai — A2 

A6 + A7 A4+A5+C Ae — A7 
y Ai — A2 Ae — A7 A1 + A2 — A3+C J 



VR3y 



(A4 - As + C) = . 



(4.2a) 



(4.2b) 



Assuming that the matrix N^j^ is invertible, we observe that R^ = 0, which implies that 
(0J02) = (020i)- This latter condition can be satisfied in two ways, if v\ 7^ 0. The first 
possibility is to have ^ = 0, which amounts to the non-breaking of the CPl symmetry by the 
vacuum. The second possibility is to have v\ = 0, with ^ 7^ and possibly t>^ 7^ 0. However, 
^ and V2 can be set to zero by an SU(2)l gauge rotation, giving rise to a CPl-invariant 
vacuum. Hence, the neutral vacuum solutions arising from an invertible matrix N^j^ do not 
break the discrete CPl symmetry and so do not lead to topological defects, such as domain 
walls. We therefore turn our attention to situations where the determinant of the matrix N^j, 
is singular, thanks to specific choices of the Lagrange multiplier (. 
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4.1.1 Neutral Vacuum Solutions from a Singular Matrix N 



In order for the matrix N^iy to have no inverse under the CPl symmetry, we require that 
the determinant of N^j, vanishes. This is guaranteed by setting the the expression in f l4.ip to 
zero. We find four possible solutions for the Lagrange multiplier, three attributed to f l4.2ap 
and one attributed to (14.2bp . However, as we have previously seen for the other symmetries 
studied so far, since the RHS of f l4.2ap is in general a non-zero vector in this case, unless 
/ij = = and Re(m^2) = 0^ this matrix equation is overdetermined. Unless the parameters 
777,^2) &iid the quartic couplings Xi^2,...,7 satisfy an unnatural fine-tuning relation, the 
matrix equation f l4.2ap becomes incompatible for the Lagrange multipliers that result from 
requiring that the matrix of f l4.2al) is singular. We therefore reject these three possible 
Lagrange multipliers and focus on the single Lagrange multiplier solution to (I4.2bp : 

C = -A45 . (4.3) 

This choice of ( lifts the constraint R^ = 0, which resulted from a non-singular matrix N^j,. 
As consequence, the CP-odd phase ^ can be non-zero in general, thus triggering spontaneous 
breakdown of the CPl symmetry after the electroweak symmetry breaking. This phenomenon 
is called spontaneous CP violation in the literature [20| 136]. 

Substituting the value of the Lagrange multiplier ( in (14. 3 p into the matrix equa- 
tion (]4.2ap . we can calculate the individual components of the 4- vector R'^. These are given 
by 

R° = — I [A5(2A2 — A345) + A67A7] fll + [A5(2Ai — A345) — A6A67] 

+ [A12A67 - (Ai2-A345)A67] "^?2} ) (4-4a) 
R' = ^ { (A345A7 - 2A2A6) /i? + (A345A6 - 2A1A7) fil + (4A1A2 - mj^ } , (4.4b) 

= ^ { [A5(2A2 + A345) - A67A7] /ii + [A6A67 - A5(2Ai + A345)] f^l 

+ [A12A67 - (A12 + A345)A67)] "^L } , (4-4c) 

with 

A = As (4A1A2 - - 2AiA? - 2X2X1 + 2A345A6A7 . (4.5) 
From f l4.4ap and (14. 4c p . we can now calculate, by means of (12. 3p . the VEVs for the bilinear 
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field expressions: 



(0^02) 



(2A2A5 — Ay) + (AeAr — A345A5) + (A345A7 — 2A2A6) 



m 



12 



A5 (4A1A2 — AI45) — 2A1A7 — 2A2AQ + 2A345A6A7 
(AeAy — A345A5) fxl + (2A1A5 — Ag) + (A345A6 — 2A1A7) 



12 



A5 (4A1A2 — AI45) — 2A1A7 — 2A2Ag + 2A345A6A7 



> 



> 



(4.6a) 



(4.6b) 



In order to fix the remaining undetermined component R^, we impose the neutral vacuum 
condition (]2.9bp on the 4-vector R^, i.e. R^ has to be a null vector. In this way, we find for 
the second component R^ that 



R^ 



4^^ 



(2A2A5 — Ay) fil + (A6A7 — A345A5) fil + (A345A7 — 2A2A6) 



m 



12 



X 



(A6A7 — A345A5) /i^ + (2A1A5 — Ag) fxl + (A345A6 — 2A1A7) 



m 



(A345A7 - 2A2Ag) fil + (A345A6 - 2A1A7) fil + (4A1A2 - Xl^^) 



12 



m 



(4.7) 



12 



1/2 



After determining all the components of R^ and comparing them with fl2.42al) and fl2.42bp . 
it is straightforward to find the vacuum manifold parameters for the Lagrange multiplier 
solution ( given in fl4.3p . with = 0. These are given by 



2 (2A2A5 — Ay) fJ'l + 2 (AgA7 — A345A5) /i2 + 2 (A345A7 — 2A2Ag) 



12 



A5 (4A1A2 — AI45) — 2A1A7 — 2A2Ag + 2A345A6A7 



2 (A6A7 — A345A5) /ii + 2 (2A1A5 — Ag) /i2 + 2 (A345A6 — 2A1A7) 



m 



12 



cos{ 



A5 (4A1A2 — AI45) — 2A1A7 — 2A2Ag + 2A345AgA7 



(4.8a) 

(4.8b) 
(4.8c) 



We note the necessary condition < |cos^| < 1, for obtaining spontaneous electroweak 
breaking of the CP symmetry in the CPl-invariant 2HDM. 

In order for the above extremal solutions to represent local minima, we require that the 
Hessian of the CPl-invariant potential be positive definite when evaluated at the extremal 
points. The Hessian with respect to V2 and ^ for the CPl-invariant potential has the 
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elements: 



Hn 

Hl2 
Hi3 



H23 
H33 



-fil + 3AiK)2 + i [A345 + 2A5 cos^ (v^^y + 3A6t^i%° cose , 



[A345 + 2A5 COS^ e] '^1^2 + COS e 



-f 2 sm 4 



2A5t;^°cose - raj, + ^AeK)^ + ^ArK)' 



H22 = -f^l + 3X2{v'^f + - [A345 + 2A5 cos' e] {v^if + SXjvy^ cose , 



2X,vy, cose - + ^AeK)' + ^A7(t;°)' 



-Vi sm ^ 
-A5K)'K)'cos2e 



-mL + ^AeKr + iA^K)^ 



f °f 2 cos e • 



(4.9a) 
(4.9b) 

(4.9c) 
(4.9d) 
(4.9e) 

(4.9f) 



It is difficult to obtain compact analytical expressions in terms of the set of potential param- 
eters {yuf , /i25 ""^125 -^15 -^25 A3, A4, A5, Ae, A7}, so the positivity of the symmetric H matrix can 
only be checked numerically for a given set of input parameters. This procedure forms part 
of our numerical analysis in Section 15.1.21 



4.1.2 CPl Topology 

The topology of the CPl-invariant 2HDM potential is very similar to the Z2-symmetric case 
discussed in Section |3l In the symmetric phase of the theory, the total symmetry group of 
the potential is Gcpi = CPl (g) SU(2)l ® U(1)y ^ Z2 x 5"^ x 5*^ Here we have used the fact 
that CPl is homeomorphic to Z2. After electroweak symmetry breaking [cf. flS.lSp ], the CPl 
symmetry breaks into the identity I, so the unbroken group is Hcpi = I ® U(l)cm — "S*^. In 
the $-space, the vacuum manifold is then given by the set 



GcPi/HcPi ^ Z2 X 53 



(4.10) 



This vacuum manifold is homeomorphic to that of the Z2 HP symmetry and we conclude 
that Ho [A^§^^] 7^ I. This implies that the CPl-invariant 2IIDM has a domain wall solution, 
which is studied in Section O 

The construction of the vacuum manifold in the $-space proceeds in a rather analogous 
manner. As stated in fl2.51al) and fl2.51bp . a general point of the vacuum manifold due a CPl 
symmetry is given by $ = A^^^^ $0) where $0 is defined in terms of the non-zero VEVs 
and the CP-odd phase e given in fl4.8al) . (14.8bp and fl4.8cl) . respectively. In the 8-dimensional 
Majorana $ space, the HP and CP transformation matrices Ai^^ of (12. 5 lap and fl2.51bl) 
have 7± = T± = cr^. Ignoring gauge transformations, there are two distinct neutral vacuum 
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solutions: 




(4.11) 

in the gauge basis, where > 0. Finally, it is worth mentioning that under the additional 
parameter restrictions = = A7 = 0, the phase ^ takes on the special value ^ = f • 
Given the freedom of reparameterization 02 — ^ i4>2 [HE], the CPl vacuum manifold coincides 
with the one of the Z2 vacuum manifold in this case. 



4.2 CP2 Symmetry 

The discrete CP2 symmetry of the 2HDM is defined as follows: 



')2 



"1 5 



'1 • 



Using the CP2 parameter restrictions of Table[8], we derive the VEVs of 0i,2 by considering the 
two conditions fl2.9al) and fl2.9bp . As before, we examine the two distinct cases: (i) det[N^j,] ^ 
and (ii) det[N^^] = 0. 

To start with, we first calculate the determinant of N^j,, which may be conveniently 
expressed as follows: 



det[N 



(2Ai + A3 - [(2Ai - A3 + C)((A4 + 0' - |Ae 

-4|A6|'(A4 - i?5 + + 8/6(/5i?6 - ^5/6 



(4.12) 



Then, for the CP2-invariant 2HDM potential, the extremization condition N^^.R'' = gives 
the two equations: 



(2Ai + A3-C)R° 
2Re(A6) \ / \ 



I A4 + Re(A5) + C -MAs) 

-Im(A5) A4 - Re(A5) + C -2Im(A6 
Y 2Re(A6) -2Im(A6) 2Ai - A3 + C / 



R2 
VR3y 







(4.13a) 



(4.13b) 
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Now, if the matrix N^,^ is invertible, the components of R'^ are found to be 

Ri = r2 = r3 = . (4.14b) 

Since only the component R° is non-zero, this result is not compatible with the neutral 
vacuum condition (12.9bp . with R'^R^ = 0, unless fif = 0. This is not a viable scenario, 
since fj'g = 0, without electroweak symmetry breaking. For this reason, we now consider the 
second possibility of a singular matrix N^j,, with det[N^i,] = 0. 



4.2.1 Neutral Vacuum Solutions from a Singular Matrix N 

We now analyze the neutral vacuum solutions, for which the determinant of N^i, vanishes due 
to a particular choice of the Lagrange multiplier (. From fl4.13ap . we see that the singular 
solution ( = 2Ai + A3 is not compatible, unless fif = 0. Therefore, we concentrate on the 
other three possible solutions obtained from requiring the vanishing of the determinant of 
the matrix on the LHS of (14.13bp . 



Employing standard methods for solving cubic equations, we obtain the three roots: 

d 6b — 2a? a / . , ^ x 

= 6 3^ 3 ' ^'-^'^^ 

^ {l+iV?,)d (l-zv/3)(3&-a^) a 

12 3d 3 ' y ■ ) 

. (l-zy3)rf , (l+«y3)(35-a2) a 

Ca — ; — , l4.ioC) 

^ 12 3d 3 ' ^ ' 

where a, 6, c and d are defined as 

a = 2Ai - A3 + 2A4 , (4.16a) 

b = 2A4(2Ai - A3) + A^ - IA5P - 4|A6p , (4.16b) 

c = (2Ai - A3)(A^ - lAsH - 4|A6|'(A4 - R5) + ShihRe - R^h) , (4.16c) 

d= (36ab - 108c - 8a^ + 12^12^3 - 3a^'^ - 54a6c + Slc^ + 12a3c) ^ ^ . (4.16d) 



Since the matrix equation fl4.13bp is under determined, we may exploit this fact to 
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express the components and in terms of R^ as 



-Re (-^4 — R5 + C) ~ hh 

3 _ 2/5-R6 - 2/6 (A4 + -R5 + C) J.1 
4/6i?6-/5(2Ai-A3 + C) 



(4.17a) 
(4.17b) 



To determine the component R^, we impose the neutral vacuum condition R'^R^ = 0. In this 
way, we obtain that 

2^^l 



Ri 



where the parameter B is given by 



±1 



B 2Ai + As - C 



(4.18) 



B 



hR^-h{\i + R^ + C) 
ReiXA — R5 + — hh 



+ 



2/5^6 - 2/6 (A4 + R5 + O 
4/6/26 -/5(2Ai- A3 + 



+ 1 



(4.19) 



We observe that there are two possible solutions for R^, and therefore for R^ and R^, 
through ( I4.17ap and (]4.17bp . The two solutions differ by a common overall sign and they 
are topologically connected via the CP2 transformation Ccp2 given in Table [3] (see also our 
discussion below in Section l4.2.2p . Considering only the positive solution of R^, the vacuum 
manifold parameters for f ^ = are calculated to be 



2/i? 



,2Ai + A3-C 



^ 1 2/5/^6 -2/6(A4 + /25 + C) 
B 4/6/?6-/5(2Ai-A3 + C) 



tan^ 



,2Ai + A3-C 
(A4 + /?5 + C)/6 - 



2/6(A4 + /?5 + C) 

B 4/6/?6-/5(2Ai-A3 + C) 



^ _ ^ 2/5/^6 



KRf 



hh — (A4 — /?5 + 0-^6 



(4.20a) 

(4.20b) 
(4.20c) 



Note that the negative solution of R^ is obtained by interchanging Vi V2 and shifting 

It is important to remark here that the phase ^ in f l4.20cp does not signal spontaneous 
breaking of the CP symmetry [37]. Within the bilinear scalar-field formalism, it is not 
difficult to see that under a unitary rotation of the Higgs doublets 01^2; which induces an 
orthogonal rotation to the 'spatial' components R^'^'^, the matrix equation fl4.13bp remains 
form invariant. In particular, one can always find an induced orthogonal rotation, such 
that the matrix on the LHS of fl4.13bl) becomes diagonal [57]. It is obvious that in this 
diagonal basis, the transformed quartic couplings Ae,? vanish and Im A5 = 0. This result is 
identical to the one found previously in [58], which is based on the construction of all possible 



42 



Quantity 


Ci 


C2 


C3 




0.372 


0.349 


0.340 




0.305 


0.261 


0.060 




-1.17 


0.343 


0.971 




-0.0578 


-0.0474 


-0.0297 



Table 10: The numerical values for the vacuum manifold parameters and potential value at 
the extremal points for the parameter set {fif, Ai, A3, A4, A5, \q} = {1, 8, 1, 3, 1 — 2i, 1 — 2i}, in 
arbitrary mass units. 



Jarlskog-like [59l EH] , Higgs-basis independent CP-odd invariants |6Tl [62] (for a recent review, 
see [53]). 

For illustration, we display in Table [10] the numerical values of the vacuum manifold 
parameters for Ci,2.3 in a CP2-invariant 2HDM, where 

{/i?,Ai,A3,A4,A5,A6} = {1,8,1,3,1-22,1-22}, (4.21) 

in arbitrary mass units. This particular set of parameters is chosen so as to satisfy the 
CP2 convexity conditions of Table [9] The values of the three Lagrange multipliers are: 
(i = —0.295, (2 = —4.09 and Cs = —16.6. In order to determine whether the three extremal 
points presented in Table [TO] are local minima, we need to analyze the positivity of the Hessian 
matrix H. 

The Hessian for the CP2-invariant 2HDM potential is a 3 x 3 symmetric matrix, with 
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the elements 
Hii 



H 



12 



Hi3 

H22 



-/i? + 3Ai(t;°)2 + i (^A34 + i?5 cos 2^ - h sin 2^) (v^^)^ 
+ (^Rq cos e - Jg sin , (4.22a) 

A34 + i?5 cos 2^ - h sin 2^) + ^ [(t;?)^ - (t;^)'] (i?6 cos ^ - h sin e) , 

(4.22b) 

(^i^gsin^ + Jecose) , (4.22c) 



H23 
H33 



-(^i?5 sin 2^ + /5Cos2ejt;5^(f^ 
-fil + 3Ai(t;°)' + i (A34 + R5 cos 2^ - h sin 2^) (t;?)' 

- (^Re cos e - Je sin , (4.22d) 

i^esin^ + /ecos^^ , (4.22e) 



- ( R5 sin 2^ + h cos 2^ ) (t;?)^^;" 



R5 COS 2^-/5 sin 2n{v^^ 



,0\2/„,0n2 



:^^2° 



.0n2 



.0n2 



^-Rg sin e — /e cos 



(4.22f) 



We can numerically check the positivity of the matrix H. In this way, we find that for a 
convex CP2- invariant potential with input parameters as given in f l4.2ip . only the Lagrange 
multiplier (i represents a local minimum, which is a global minimum. As we will see below, 
this global minimum has a twofold degeneracy, as a consequence of the CP2 symmetry. 

4.2.2 CP2 Topology 

In the symmetric phase of the theory, the total symmetry group of the CP2-invariant 2HDM 
potential is Gcp2 = CP2 (g) SU(2)l ® U(1)y = Z2 ® Ha ® SU(2)l ® U(1)y, where U2 is the 
permutation symmetry 0i 4-7- 02- To be specific, we have used here the isomorphism [23]: 
CP2 = 22 0112, which is evident in the Z2-constrained Higgs basis where Ag = A7 = 15 = 
0. After electroweak symmetry breaking [cf. fl3.18p ]. the permutation symmetry 112 remains 
intact, so the residual unbroken group of CP2 is Hcp2 = ^2 ® U(l)em- As a consequence, the 
vacuum manifold A^§^^ in the <I>-space has the topology of the coset space: 



A^$^^ — GcP2/HcP2 



(4.23) 



The vacuum manifold A^^^^ is homeomorphic to that of the Z2 HF symmetry, thus having 
a non-trivial zeroth homotopy group Ho [A^^^^] = IIq [Z2] 7^ I. This implies that the CP2- 
invariant 2HDM has a domain wall solution, which we analyze in Section O 
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An arbitrary point $ of the vacuum manifold due to a CP2 symmetry may be obtained 
with the help of fl2.51ap and (]2.51b|) . i.e. $ = A^^^^ $0) where $o is defined in terms of the 



non-zero VEVs 2 given in f l4.20ap , fl4.20bp and ^ in fl4.20cp . The HF and CP transformation 
matrices Ai'^'^ of (12.5 lap and (12.51bp are 7+ = T_|_ = and 71 = T_ = icr^, respectively. 
From the action of these transformation matrices on $0) "we find that the vacuum manifold 
is comprised of two disconnected sets. The elements within each set are related by gauge 
rotations U±. Two representative vacuum manifold points from each set are 









(4.24) 



and 





(4.25) 



where we used the freedom of the gauge rotations U-t, in order to adjust the neutral compo- 
nent of 01 to be positive. 

4.3 CPS Symmetry 

The CPS symmetry is a continuous CP symmetry and is defined by the transformations 

01 — > (p'l = COS (/){ + sin 6 02 , 

02 — ^ 02 = ~ sin + cos 9 02 , 

where 9 G [0, vr). As before, we first consider the case det[N^i,] 7^ 0. The determinant of N^^, 
is given by 

det[N^,] = (2Ai + A3 + C) (2Ai - A3 - C)' (2A4 + A3 - 2Ai - C) • (4.26) 

For the CP3-invariant 2HDM potential, the extremization condition N^j.R'^ = leads to 
four separate equations: 

(2Ai + A3-C)R° = 2/i2, (4.27a) 

(2Ai-A3 + C)R' = 0, (4.27b) 

(2A4 - 2Ai + A3 + C)R' = 0, (4.27c) 

(2Ai - A3 + C)R^ = . (4.27d) 
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Based on the assumption that N^,^ can be inverted, we find that all "spatial" components 
j^i,2,3 vanish. Like in the CP2 case, the condition for having a neutral vacuum restricts the 
remaining component R° to be zero as well, which can be naturally fulfilled, only if /x^ = 0. 
In such a scenario, one has = ^ and so absence of electroweak symmetry breaking. 
Therefore, we now investigate the case where det[N^y] = 0. 



4.3.1 Neutral Vacuum Solutions from a Singular Matrix N 

From the system of equations fl4.27al) - fl4.27dp . it is easy to see that there are only two 
compatible singular solutions of N^,y for the Lagrange multipliers: 

Ci = -2A1 + A3, (4.28a) 

C2 = 2Ai - A3 - 2A4 . (4.28b) 



Let us first consider the solution (i = — 2Ai + A3. In this case, only the components R° 
and R^ of the 4- vector R^ are determined as 

R° = ^ , (4.29a) 
R2 = . (4.29b) 

Instead, R^ and R^ are free parameters, which are constrained by the neutral vacuum condi- 
tion: R'^R^ = 0. Specifically, the latter condition gives rise to the constraint: 

i^'r + = [0 ■ (4.30) 

The constraint R^ = implies that ^ = rnr, where n is an integer. In terms of the vacuum 
manifold parameters fj'g; we have the general solution 

= sin^, = cos^ , (4.31) 

vAi VAi 

where ^ = tt and 6 G [0, tt). The free angle 6 is associated with a massless 'CP-even' Gold- 
stone boson, as can be verified independently from the analytical results presented in [TB] . 

In order for the extremal point given in fl4.3ip to be a local minimum, we require that 
the elements of the Hessian matrix H for the CP3-invariant 2HDM potential be positive. The 
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elements of the symmetric matrix H read: 

Hn = -fxl + + ^ [A34 + (2Ai - A34) cos2e] (t;2°)' , (4.32a) 

H12 = [A34 + (2Ai - A34) cos 2^] , (4.32b) 

Hi3 = - (2Ai - A34) v',{v^^y sin 2^ , (4.32c) 

H22 = -/i? + 3Ai(t;°)2 + ^ [A34 + (2Ai - A34) cos2e] {v'.Y , (4.32d) 

H23 = - (2Ai - A34) KM sin 2e , (4.32e) 

H33 = - (2Ai - A34) (t^?)'(t^2)' cos 2e . (4.32f ) 

Then, the conditions for positivity of H are simply given by 

fil > , A34 > 2Ai . (4.33) 

Note that the second condition in fl4.33p is supplementary to the two conditions given in 
Table [9] for ensuring a convex CP3-invariant 2HDM potential. This local minimum has the 
potential value 

Vo ^ . (4.34) 

Let us now investigate the second singular solution (2 = 2Ai — A3 — 2A4. In this case, 
we obtain 

rO = , , (4.35a) 

2(A3 + A4) ' ^ ' 

Ri = r3 = 0. (4.35b) 

The component is constrained by the neutral vacuum condition fl2.9bp imposed on R'^, 
i.e. R'^R^ = 0, from which we find that 

R^ = ± R° . (4.36) 



Taking the constraints f l4.35ap . f l4.35bl) and f l4.36p into account, we derive the vacuum man- 
ifold parameters 

with v' = /U1/A/A34. The conditions for this neutral vacuum solution to be a local minimum 
result from the positivity of the Hessian matrix H given in (I4.32ap -( l4.32fp . These conditions 
are 

fil > , 2Ai > A34 . (4.38) 
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These conditions are, in general, not guaranteed solely by the convexity conditions for the 
CP3-invariant potential stated in Table M Direct comparison of these minima conditions 
with those for the Ci solution in fl4.33p shows that both local minima cannot co-exist. It 
depends on the relative values of 2Ai and A34 which solution becomes the local minimum, 
and this will then be the global minimum as well. The value of the potential arising from 
the second solution (2 is easily evaluated to be 

Vo = - ^ . (4.39) 
ZA34 

In the following, we will analyze the topology resulting from the two neutral vacuum solutions 
given in fl4.3ip and fl4.37p . respectively. 

4.3.2 CPS Topology 

It is interesting to discuss the vacuum topology of the CP3-symmetric 2IIDM for the two so- 



lutions obtained by means of the Lagrange multipliers Ci and (2, given in (I4.28ap and (I4.28bp . 
respectively. 

We first note that the total symmetry group of the CP3-symmetric 2IIDM potential 
is GcP3 = CP3 O SU(2)l O U(1)y ~ Z2 x 5^ x 5^ x S\ since CP3 ^ CPl ® S0(2). This 
means that the CP3 group is equivalent to the combined, as well as independent action of a 
standard CPl transformation and a S0(2) HF rotation in the {(pi , ^2) field space. 

Let us now consider the neutral vacuum solution obtained by the Lagrange multiplier (i. 
After electroweak symmetry breaking, the total symmetry group Gcps breaks into the residual 
group Hppg = CPl ® I ® U(l)em- This can be easily seen, since ^ = in this scenario and 
so the CPl symmetry remains intact, whereas the SO (2) HF symmetry gets spontaneously 
broken to the identity I. As a consequence, the vacuum manifold in the $-space is determined 
by the coset space 

-Mf' = GcP3/Hg^3 - S'xS\ (4.40) 

Since Hi [M2^^] = Hi [S^] ^ I, we conclude that the CP3-invariant 2HDM related to the 
Lagrange multiplier (i has a vortex solution which is analyzed in detail in Section [51 Using 
the result of (14.310 . the transitive action of the transformation matrices of (I2.51ap and (]2.51bp 
result in the general points on the vacuum manifold: 

^1 = ^1 ° I , '^2 = ^1 ° I , (4.41) 

V2 \^ ycos9 J V2 \^ (-l)"i;sin^ J 

where v = Hi/y/Xl and 6 E [0,7r). There is a relative minus sign for odd n, but this can be 
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absorbed by redefining 9 as tt — 6. 

We may now determine the vacuum manifold of the CP3-symmetric 2HDM associated 
with the second Lagrange multipher solution (2 fl4.28b|) . In this case, the total symmetry 



(2) 

group follows a different breaking pattern and the little group is H^pg = CPl®SO(2)(8)U(l)c 
i.e. neither of the two symmetries CPl and SO (2) are broken. In order to see this, we may 
consider an SO (2) rotation of the vacuum manifold point given in fl4.37p . yielding 









(4.42) 



The phase 6 can always be removed by a U(1)y hypercharge rotation, which is a manifestation 
of the fact that the SO (2) symmetry is not broken, after electroweak symmetry breaking. 
Moreover, one could reparameterize the second Higgs doublet 02 as ±202, in order to render 
both VEVs of 01^2 real. Since such a reparameterization does not induce any additional 
phase in the real quartic couplings of the CP3-invariant 2HDM potential, we conclude that 
the CPl symmetry is not broken as well. Thus, the vacuum manifold determined by the 
coset space Gcps/H^^pg is homeomorphic to S^, exactly as in the SM. Consequently, there 
are no non-trivial topological defects in the 2HDM scenario related to the second Lagrange 
multiplier solution (2- 



5 Topological Defects in the 2HDM 

Using our analysis of the six accidental symmetries of the 2HDM conducted in Sections E] 
andlU we will now study the topological defects associated with the spontaneous symmetry 
breaking of each of the six accidental symmetries studied. As shown in Table [TTl we find 
that there are three domain wall, two vortex and one global monopole solutions due to the 
additional symmetries of the 2HDM, possibly posing significant cosmological implications 
for the model. A comprehensive introduction to the properties and formation of topological 
defects is given in [?T]. 

In our study of the topological defects, we assume that the VEVs of the two Higgs 
doublets 01.2 are still assigned at and after the electroweak symmetry breaking, such that 
{v^f + {v^f = vl^, where vsm ~ 246 GeV is the VEV of the SM Higgs doublet. Due 
to the complexity of the differential equations that result from the 2HDM Lagrangian for 
each symmetry, our study of the scalar functions involved is carried out numerically using 
gradient flow techniques which involve minimizing the energy of a configuration on a finite 
grid with initial conditions that have the appropriate boundary conditions. This is done by 
defining an energy functional E = E(/i, . . . , /„), where /i, . . . , /„ are the functions defining 
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U(1)PQ 
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I 




Vortex 


S0(3)hf 


S0(3)hf 


S0(2)hf 




Global Monopole 


CPl 


CPl ~ Z2 


I 


Z2 


Domain Wall 


CP2 


Z2®n2 


n2 


Z2 


Domain Wall 


CPS 


CPl ® S0(2) 


CPl 




Vortex 



Table 11: Breaking patterns of the total symmetry group Ghf/cp in-to the little group Hhf/cp; 
after the electroweak symmetry breaking SU(2)l (S) U(1)y — ?■ U(l)j,jjj. The fourth and fifth 
columns show the topology of the vacuum manifold A^^^^^^ and the associated topological 
defect, for each of the six accidental HF/ CP symmetries of the 2HDM. 

the topological solutions, and then by solving the first order diffusion equation fk = —jf' 
k = 1, . . . , n. 

5.1 Domain Walls 

We begin our discussion of topological defects with domain walls, which have long been 
known to have severe consequences for the evolution of the Universe should they form at a 
symmetry breaking phase transition in the early Universe, since they can come to dominate 
the Universe's energy density [M]. Various mechanisms to reconcile this undesirable nature 
of domain walls with current observations have been discussed, such as the restoration of 
the broken discrete symmetry and subsequent evaporation of the domain walls at a later 
phase transition the use of a period of exponential inflation to dilute the concentration 
of domain walls [BS] and the symmetry of the model being only an approximate discrete, 
exponentially suppressing domain wall energy density [HZl ESI EH] • 

The present study of domain walls does not attempt to analyse the cosmological impli- 
cations, which will be presented in a future publication, rather it focuses on presenting an 
overview of the typical domain wall solutions and analysing whether or not the energy per 
unit area of the domain wall can be made to be vanishingly small for specific valid parameter 
choices. 

5.1.1 Z2 Domain Walls 

From our analysis in Section I3.H the 2HDM potential that is invariant under the HF Z2 
symmetry can exhibit a disconnected vacuum manifold, the components of which are not 
linked by the gauge symmetries of the theory, provided both VEVs of the Higgs doublets 
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(j)i^2 that create the neutral vacuum global minimum solution are non-zero, i.e. 7^ and 
V2 = 0- This scenario is only apparent within the Z2 invariant 2HDM when the determinant 
of the matrix N^,^ vanishes, as shown in Section 13.1.11 

Let us now analyze an one- dimensional, time-independent kink solution for the Z2 
symmetry. In order to find such a solution, we will use an ansatz for the two Higgs doublets 
given by: 

0i,2(x) = i= f ° I , (5.1) 

where the coordinate x describes the spatial dimension perpendicular to the plane of the 
domain wall. Using this ansatz, the energy per unit area of the system is 



E = dx £(01,02) , (5.2) 
^ — 00 

where the energy density for the general 2HDM is given by: 

£:(0i,02) = (V0l) ■ (V0i) + (V0^) ■ (V02) + V(0i,02) + Vo , (5.3) 

where V is the 3-dimensional gradient operator, expressed in the relevant coordinate system. 
Moreover, Vq is introduced to normalize the potential contribution to the energy density to 
have a zero value at the global minimum. The energy density for the Z2 invariant 2HDM is 
given by 



dx J 2 \dx J 2^' ' 2' 
+ hiv'A^r + \^2v'2{xr + i (A34 - IA5I) v',{xfv',{xY + Vo . (5.4) 

To simplify our study, we introduce the dimensionless quantities 

Vi2{x) . A2E 

X = /i2X , ■01,2(3^) = — ■ ; E = — 3- , (5.5) 

in order to rescale the energy per unit area of (15. 2p to be dimensionless. Here, we introduce the 
convention that " represents a dimensionless quantity. Performing these rescalings leaves the 
dimensionless Z2 energy density, denoted correspondingly as dependent on the following 
parameters: 

2 lA \ A34— IA5I (12 /p. „N 



H2 A2 2A2 VA2 

Also, the vacuum manifold parameters of fl3.13ap and (I3.13bp . which are the boundary 
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conditions on the fields f°2(a;), are rescaled, such that 



hm v^{x) 



hm ■02(3;) 




(5.7) 



The equations of motion for the two rescaled fields Vi2i^) found to be 



dHl 
dx"^ 



(5.8a) 
(5.8b) 



As no analytical solutions are known for this particular system of ordinary differential equa- 
tions, we proceed by gradient flow techniques to minimize the energy per unit area. To make 
this possible, we truncate the interval of integration of (15. 2 p from (—00,00) to [— -R, i?], en- 
suring that R is chosen to be much larger than the width of the kink. By making the range of 
integration symmetric about x = 0, we break the translational symmetry usually exhibited 
by kink solutions. 

In order to perform the numerical analysis, a particular parameter set {/i^. A, (7} must 
be chosen that satisfies the constraints for a bounded-from-below global minimum, as these 
are given in f l3.15p and Table O For convenience, we state these results in terms of the 
rescaled parameter set: 



\ 2 2 A 

A>fi', g < II < - , 

9 



9 > 



(5.9) 



Additionally, in order to satisfy that the sum of the squares of the two VEVs is "^gj^, 
require that the VEV scale factor 77 have the value given by 



we 



V2 



A 



9' 



VSM ■ 

9 



(5.10) 



Here, we make the observation that a value of rj can always be found that ensures condition 
f l5.10p is met for any parameter set {fi^, X,g}, provided that the members of the parameter 
set remain non-zero, finite and satisfy conditions (15. 9p . 

We present two typical solutions in Figure[2]for the parameter sets {/i^. A, g} = {1, 1, 0.5} 
and {1.5, 1,0.5}. We also show the general form of the dimensionless energy E in Figure [3] 
and directly compare several different solutions in Figure IH as a function of /i^. From Figures 
[3] and m we see that as /i^ approaches its lower bound, /i^ — )■ g, the dimensionless energy 
approaches a finite value. In the limit g 0, we find that this finite value is the familiar 
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Figure 2: Plots of Vi{x), V2ix) and the dimensionless energy density £{x) for two different 
valid parameter sets of the Z2 invariant potential. The parameter sets used are {1, 1,0.5} 
(LHS) and {1.5, 1,0.5} (RHS), and the region of integration has R = 15. 



1.0 




^2 



Figure 3: Numerical evaluation of the dependence of the dimensionless energy E as a function 
of n^, for the Z2 invariant potential. Here, we use the fiducial values A = 1 and g = 0.5. 
Convexity of the potential and global minima conditions for these values require that fi"^ G 
(0.5,2.0). 



value |"\/2 [cf. (1.20) in |12]], and the kink width decreases and becomes smalL Conversely, 
we see that as /i^ approaches its upper bound, i.e. as /x^ — )■ ^, the dimensionless energy E 
tends towards zero, the kink width increases and the energy density becomes delocalized. 
Therefore, the dimensionless energy can be made vanishingly small for appropriate choices of 
the parameter set. This is a feature that can be exploited to avoid domain wall domination 
by making the mass per unit area of the walls ultra- light. 
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Figure 4: Plots comparing various Vi{x) (LHS plot) and V2{x) (RHS plot) curves for the Z2 
invariant potential by fixing A = 1 and g = 0.5, and allowing /i^ to vary, and the region of 
integration R = 15. 



5.1.2 CPl Domain Walls 

From our analysis in Section 14.11 the 2HDM potential which is invariant under the CPl 
symmetry can exhibit a disconnected vacuum manifold, the components of which are not 
linked by the gauge symmetries of the theory, provided both VEVs of the Higgs doublets 01^2 
and the relative phase between the doublets that create the neutral vacuum global minimum 
solution are non-zero, i.e. 7^ 0, = and ^ 7^ 0. This spontaneous violation of CP 
is only apparent within the CP 1-invariant 2HDM when the determinant of the matrix N^,y 
vanishes, as shown in Section 14.1.11 However, as we have shown in Section 14. H a neutral 
vacuum solution with a CP-conserving global minimum is also possible, where det[N^j,] 7^ 0, 
i.e. a global minimum with 7^ 0, = and ^ = 0. Therefore, during our numerical 
analysis we are careful to choose parameter sets that give global minimum neutral vacuum 
solutions with spontaneous CP violation. 

Let us now investigate an one-dimensional, time- independent kink solution for the CPl 
symmetry. In order to find such a solution, we will use an ansatz for the two Higgs doublets 
given by 

1 / \ 1 / \ 

where the coordinate x describes the spatial dimension perpendicular to the plane of the 
domain wall. The energy per unit area associated with the kink solution is again given by 
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fl5.2p . where the energy density for the CPl-invariant 2HDM is given by 



2\dx J 2\dx J 2 ' \dx J 2"^' ' 2 
+ 1Ai<(x)^ + \x2v',{xY + \ (A34 + As cos 2e(x))t;?(x)\°(x)2 

+ ( - + IKv^x)' + hy,{xfy^,{xy,{x) cose(x) + Vo . (5.12) 

By rescahng (15 ■2p to be dimensionless for the CPl energy density, we again make use of the 
dimensionless quantities of (15. 5p . Performing these rescahngs leaves the dimensionless CPl 
energy density dependent on the following parameters: 

2/^1 2 ^12 X Ai A34 fr , r « '7^ 

;x = — , m = , A = — , 5(34 = — , gk = Y- fc = 5, 6, 7j , i] - 



jd,2 A2 A2 A2 '\/A2 

(5.13) 

It is also useful to introduce the parameter g = (734 — g^. The parameter set for the CPl- 
invariant model then reduces and becomes {fi'^,m'^, X, g^^, g^, ge, gj}. Similarly, the vacuum 
manifold parameters f°2 and ^ of (I4.8ap . f l4.8bp and (I4.8cp . which are the boundary conditions 
on the fields v^2i^) ^i^)^ rescaled, so as to give 



--±- ' ^ ^ V (4A - ^2) _ 2Xgl - 2(7? + 2gg,gr ^ ' 

y .o.^^ _ / 2 (2A(75 - gl) + 2 ((76^7 - m) /^^ + 2 (^(76 - 2A(77) m2 

^iim^.2(x) - y (4A - - 2A(7| - 2(7? + 2^(76^77 ' ^ ^ 

hm m = iarccos f ' ^^^'"^^ ^^^"^^^ ) . (5.14c) 
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Figure 5: Plots of v\{x), V2{x), ^{x) and the dimensionless energy density S{x) for two 
different valid parameter sets of the CPl-invariant potential. The parameter sets used are 
{1,0.1,1,2.5,1,0,0} (LHS) and {1,0.1,1,2.5,1,-0.15,0.15} (RHS). The region of integra- 
tion has R = 15. 



The equations of motion for the three rescaled fields vl2{^) and ^(x) are 



dx"^ 



dx'^ 



dx'^ 



+ A({;°)2 + ^ ((734 + COS 20 {vlf + ^g^vy^ cos ^ 



5? 



- '-gMf 



V2 cos ^ . 



-1 + + 1 ((734 + (75 COS 20 (i)?)' + '^gjv'ivl cose 



3 

2' 



.ocose + (|)%° 



^\ fd^ 

dx J \ dx 



V1V2 sin^ 



(75i)i%2° cose - + ^^76(^)?)' + lg7{v^2? 



(5.15a) 



(5.15b) 



(5.15c) 



As with the Z2 domain wall study, no analytical solution is found to these equations of motion 
and we therefore proceed by gradient flow techniques. Due to the number of individual 
parameters one may tune within the confines of the CPl convexity and minima conditions, 
relationships between the parameters are in general complicated and so we end our CPl 
domain wall study by presenting two typical solutions in Figure [51 However, we do note 

TT 

two cases determined by specific choices of the parameter set. For the case lim ^{x) = — , 

£—^■±00 2 

which is guaranteed if geiviY + g7{v2)'^ = 2m^, the CPl symmetry domain wall reverts back 
to a Z2 style domain wall by use of the reparameterization 02 i4>2, as discussed in Section 
14.1.21 An explicit example can be seen for the Z2 symmetry potential parameters constraints. 



m 



12 



A, 



A7 = 0. Also, the dimensionless energy E can be made vanishingly small for 
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certain valid choices of the parameter set that leave the limit of f l5.14cp still finite but ^ 1. 
This is consistent with the case lim ^(x) = in which spontaneous CP violation ceases and 
subsequently there is no domain wall solution, as discussed in Section I4.1.2I 



5.1.3 CP2 Domain Walls 

From our analysis in Section \4.2\ the 2HDM potential which is invariant under the CP2 
symmetry can exhibit a disconnected vacuum manifold, the components of which are not 
linked by the gauge symmetries of the theory, provided both VEVs of the Higgs doublets 0i 2 
that create the neutral vacuum global minimum solution are non-zero, i.e. 7^ 0, = 0. 
This scenario is only apparent within the CP2-invariant 2HDM when the determinant of the 
matrix N^j, vanishes, as shown in Section 14.2.11 

We now investigate an one-dimensional, time-independent kink solution for the CP2 
symmetry. In order to find such a solution, we will use an ansatz for the two Higgs doublets 
given by flS.lip . The energy per unit area associated with the kink solution is again given by 
fl5.2p . The energy density for the CP2-invariant 2HDM is given by 

+ \Xi{vUxr + v'i^Y) + ^(A34 + i?5Cos2e(x) - /5sin2e(a;))<(x)\°(x)2 
+ ^v^i{x)v^{x){v'^^{xf -v^{xf){RQCOs^{x) - lesin^ix)) + Vq . (5.16) 

Again, it is useful to rescale (15. 2p to be dimensionless for the CP2 energy density, and so we 
introduce 

^12(3^) - AiE 

X = fiix , Vinix) = — '■ , E = — ^ . (5-17) 

Performing these rescalings leaves the dimensionless CP2 energy density S dependent on the 
following parameters: 

gk = (for /c = 3,4,5,6) , hk = ^ (for A; = 5, 6) , i] = -7^. (5.18) 



Therefore, the parameter set for the CP2-invariant model becomes {gs, g^, g^, g^, h^, hg}. Also, 
the vacuum manifold parameters w°2 and ^ of (l4.2Uap . (I4.20bp and (l4.2Ucp . which are the 
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boundary conditions on the fields fi2(^) ^{x), are rescaled as follows: 



\ \2 + g3-C 



l_ }^ "^hge - 2he{g4 + 95 + 
B 4^6-/^5(2 -(73 + 



, as a; — )■ — oo 



Vi[X) —7- < 



^ _^ 1_ 2/155(6 - 2/16(54 + 5-5 + 



\ V2 + 53 - C/ V ^ 4/^656 - /i5(2 -53 + 



as X — )■ +00 



(5.19a) 



\ V2 + 53-C 



^ 2/^556 - 2/16(54 + 55 + 

B 4/1656 -/^5(2-53 + 



^V2 + 53-C 



^ _ }_ 2/1556-2/15(54 + 55 + 

B 4/^656 - /i5(2 - 53 + 



, as X — 7- —00 

as X — !■ +00 

(5.19b) 



lim ^(x) 



arctan 



{94 + 95 + OK - hge 
h5hfi - (54 - 55 + 056 



(5.19c) 



where the parameter B is defined as 



B 



( h596 - he{94 + 95 + C) \ ^ ( 2 /?-556 - 2/^6(54 + 55 + 



\ V5'6(54 - 55 + - hhe J \ Ah^g^ - /i5(2 - 53 + 



+ 1 . 



(5.20) 



These boundary conditions depend on the non-trivial Lagrange multiplier implemented to 
produce the neutral vacuum solution. This Lagrange multiplier satisfies the cubic equation: 

C^' + (2 - 53 + 254)C2 + (254(2 - 53) + 54 - 55 -hi- - ^hl)C 

+ {2-93){9l-9l-hl)-A{9l + hl){9,-9,) + 8K{h9,-9,h,) = 0. (5.21) 



In order to find a valid parameter set, we start by choosing parameter values that satisfy the 
CP2 convexity conditions (shown in Table [9]) and then solve fl5.2ip to find the three possible 
values of We then find the rescaled vacuum manifold parameters which correspond to each 
( solution, and determine if these solutions correspond to local minima, i.e. we require that 
the CP2 Hessian matrix H in fl4.22ap -( l4.22fp be positive definite. If they do indeed relate to 
minima, we calculate the value of the potential at these extremal points to determine which 
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Figure 6: Numerical estimates of v\{x), V2{x), ^{x) and the dimensionless energy density 
£{x), for a valid parameter set of the CP2-invariant potential. The input parameter set is 
{0.125, 0.375, 0.125, 0.125, -0.25, -0.25}. The region of integration has R = 15. 



( solution generates the global minimum. 

The equations of motion for the three rescaled fields "012(2;) and ^(x) are: 



dHl 
dx'^ 



-1 + (^?)' + 2 (9u + 9^ cos 2^ - /i5 sin 2^ {v''^? 
+ \vl{Z[v^^f - ({)°)')((76cose-/i6sin0 , 

-1 + {v^^f + ^ ((734 + 9, cos 2^ - h, sin 2^ {v^f + 
+ \vl{{vlf-m?){9^^o^i-h,sm0 , 

+ ((t)?)2-({)°)2)(^6sine+/i6cos0] . 



(5.22a) 



(5.22b) 



(5.22c) 



We can obtain numerical solutions to these equations of motion, by making use of gradient 
flow techniques. Since there is a large number of individual parameters that could vary 
within the confines of the CP2 convexity and minima conditions, relationships between the 
parameters are, in general, complicated. Instead, we present a typical example of CP2 domain 
walls, as shown in Figure [6l As with the previous two case, we observe that the dimensionless 
energy E can be made vanishingly small for certain valid choices of the parameter set, such 
as allowing g^^ and /is^e to tend to zero. 
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5.2 Vortices 



We now turn our attention to other topological defects which may form in the 2HDM, such 
as vortices. Whilst vortex solutions have been discussed in the 2HDM [70j, these vortices 
were generated by the SM gauge group, whereas the vortices we discuss here are generated 
solely by the spontaneous breaking of the U(l)pQ and CPS accidental symmetries which the 
2HDM can exhibit for specific constraints on the parameters of the potential. 

Vortices are often regarded as the most favourable topological defect, since their energy 
density does not grow relative to the background and so for sufficiently small initial energy 
densities, vortices behave benignly and can comply with current cosmological observations. 
Due to the axially symmetric, one- dimensional and typically high mass density characteristics 
of a cosmic string, strings can act as a gravitational lens [TH [72] and searches are already 
under way to detect possible vortices which may be within the current horizon, e.g. using 
precision cosmic microwave background data from experiments such as WMAP [731 17^ - 

As with our domain wall study, we do not study the cosmological implications of the 
2HDM's vortices, which we reserve for a future study, but focus on presenting an overview 
of typical solutions and determining whether or not the energy per unit length of the vortex 
can be made to be vanishingly small for specific and valid choices of the model parameters. 



5.2.1 U(l)pQ Vortices 

From our analysis in Section 13. 2[ the 2HDM potential which is invariant under the global 
Peccei-Quinn U(l) symmetry can exhibit a non-simply connected vacuum manifold provided 
both VEVs of the Higgs doublets 0i,2 that create the neutral vacuum global minimum solution 
are non-zero, i.e. 7^ and = 0. This scenario is only apparent within the U(l)pQ 
invariant 2HDM when the determinant of the matrix N^^^, vanishes, as shown in Section 13.2.11 

In order to find a two-dimensional time-independent vortex solution for the U(1)pq 
symmetry, we use the ansatz for the two Higgs doublets: 

Mr) = ^ I ° , Mr,x) = 4^ ( , ° . h (5-23) 

where the coordinate r describes the space radially outward from the core of the vortex, and 
X is an azimuthal angle which accounts for the winding of the vortex, with winding number 
n. Using this ansatz, the energy per unit length of the system is then: 

POO 

E = 27r / rdr ^(0i,02) , (5.24) 
Jo 
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where the energy density for the U(l)pQ invariant 2HDM is given by 

1 f dv^Y 1 fdvlV 2 1 2 0/ ^2 1 



+jAii.f(r)-' + hWKr)" + iA34«?(r)X(r)^ + V„ . (5.25) 

For this type of energy density, the integral of (15.241) is logarithmically divergent in r, so we 
truncate the region of integration from [0, oo) to [0, -R], where i? is a cut off radius To 
once again simplify our study, we rescale the energy per unit length of the vortex in (15.241) 
to be dimensionless by introducing the dimensionless quantities 

-0 / \ '^l,2(^) - E 

r = /i2r , v^ ^{r) = , E = - — - . (5.26) 

7] 27r?7^ 

With the help of these rescalings, the dimensionless U(1)pq energy density derived from 
(I5.24P depends now on the following parameters: 



= — , ^ = IT ^ 3 = 7^, V = ■ (5.27) 

'^2 ^^2 V^2 



Having rescaled the vacuum manifold parameters 2(^)1 we require that these approach 
their corresponding VEVs given in (I3.26al) and (13.26bp . as r — )• 00. To be precise, we impose 
the boundary conditions: 



^ =0, lim{)?(f) = W^-^, (5.28a) 



f=0 



r— >oo 





-9 




9 ' 

-r 


A- 


-^^^ 




-9' 



lim{;^(f) = 0, lim{;^(f) = \ l ^ . (5.28b) 

These conditions force "W^l^) be regular for all values of f and require ^^'(f) to be continuous 
and radially symmetric. The equations of motion for the two rescaled fields 'Oi2(^) are found 
to be 

^ = f,0(_^VA(S°)'+9ftV), (5.29a) 

^ + ?§ = *° ( - ^ + ?J + + ■ <^-'^''' 

As is typical of vortex studies, no analytical solutions to the equations of motion are found 
and so we make use of gradient flow numerical techniques. 

For our numerical analysis, we choose a particular parameter set {/x^. A, g, n} that sat- 
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Figure 7: Plots of Vi{f), V2{r) and the dimensionless energy density £{f) for two differ- 
ent valid parameter sets of the U(l)pQ invariant potential. The parameter sets used are 
{1,1,0.5,1} (^Li/Sy anc? {1.5, 1, 0.5, 1} (RHS). The cut off radius used for both plots is R = lb . 



isfies the constraints for a bounded-from-below global minimum, which are of exactly the 
same form as for the Z2 symmetry in f l5.9p . In order to satisfy that {v\Y + (^2)^ — ^sm' "^^ 
require that the VEV scale factor rj have the value given by (I5.10p . We note that a value of rj 
can always be found that ensures condition f l5.10p is met for any parameter set {/i^. A, (7,n}, 
provided that the members of the parameter set remain non-zero and finite, and satisfy the 
conditions in fl5.9p . 

We conclude our U(l)pQ vortex study by presenting two typical solutions in Figure [3 
We show the general form of the dimensionless energy E in Figure [H as a function of /i^, 
noting that the dimensionless energy tends to zero as yU^ approaches its upper limit, i.e. 
^2 _s, A We also directly compare several different solutions in Figures |9] and [TO] by varying 
/i^ and the winding number n respectively. From Figures [9] and [10] in particular, we see that 
as /i^ increases, the width of the vortex core increases and similarly, as the winding number 
increases, so does the width of the vortex core. 

5.2.2 CPS Vortices 

As discussed in Section I4.3[ the CP3-symmetric 2HDM potential can exhibit a non-simply 
connected vacuum manifold, provided a neutral vacuum global minimum solution exists 
where the sum of the squares of the two VEVs of the doublets 0i^2 is non-zero. Such a 
scenario can be realized within the CP3-invariant 2HDM, if the matrix N^j, happens to be 
singular. However, as shown in Section [1.3.1[ there are two possible neutral vacuum solutions 
that could form the global minimum solution, depending on the relative magnitudes of the 
quantities 2Ai and A34. If 2Ai > A34, we find that any possible vortex solution can be removed 
by gauge transformations, whereas for cases with A34 > 2Ai, no such gauge transformations 
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Figure 8: Numerical evaluation of the dimensionless energy E, as a function of for the 
U(l)pQ invariant potential. Here, we use the fiducial values X = 1, g = 0.5 and n = 1. 
Convexity of the potential and global minima conditions for these values require that iJ? G 
(0.5,2.0). 




Figure 9: Plots comparing various Vi{f) (LHS plot) and V2{r) (RHS plot) curves for the 
U(l)pQ invariant potential. Here, we fix X = 1, g = 0.5 and n = 1, and allow /x^ to vary. 
The cut off radius used for both plots is R = 15. 



are possible, allowing a vortex solution. Hence, we study cases of the latter type to ensure 
vortex formation in the CP3-invariant potential. 

In order to obtain a time-independent vortex solution for the CPS symmetry, the fol- 
lowing ansatz for the two Higgs doublets is used: 

Mr,x) = ° I , Mr,x) = ° I , (5.30) 

y/ Z \ v{r)cos{nx) I \ — f (r) sin(nx) 

where the coordinate r describes the space radially outward from the core of the vortex, and 
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Figure 10: Plots comparing various Vi{f) (LHS plot) and V2{r) (RHS plot) curves for the 
U(l)pQ invariant potential for various winding numbers n. Here, we fix X = 1, g = 0.5 and 
IJ^ = 1, and allow n to vary. The cut off radius used for both plots is R = 50 



X is an azimuthal angle which accounts for the winding of the vortex, with winding number 
n. Using this ansatz, the energy per unit length of the system is given by fl5.24p where the 
energy density for the CP3-invariant 2HDM is given by 

^(0 = + ^Mrf - If^lvir)' + ]xMry + Vo . (5.31) 



2 \dr J 2r2 ' ' 2"^' ' ' 4 

As before, fl5.24p is logarithmically divergent for the CPS energy density and so we truncate 
fl5.24p to a cut off radius r = R. Our study can be simplified, if we rescale the energy per unit 
length of the vortex in fl5.24p to be dimensionless and introduce the dimensionless quantities 

v(r) - E 

f = /iir , vir) = ^ , E = — ^ , (5.32) 

with r] = -y=. We then require that (f?)^ + (^2)^ = and also rj = v^u- Under the 



/AT' 

above parameter re-definitions and provided the winding number n is non-zero, the boundary 
conditions on the vacuum field v{f), which follows from fl4.29ap . become 

limv{r) = , lim v{f) = 1 . (5.33) 

These conditions force v{f) to be regular for all values of f and ensure that the dimensionful 
field v{r) approaches its VEV in the limit r — ?■ 00. The equation of motion for the rescaled 
field v{f) is 

-pi + -T^ = ^ -1 + - + ^ • (5-34) 
ctr^ r ar V / 

The above differential equation is solved numerically, by means of gradient flow methods. 
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Figure 11: Plots of v{r) (LHS plot) and the dimensionless energy density £{r) (RHS plot) 
for the CP3-invariant potential. The winding number n is varied from 1 to 5 and the cut off 
radius for both plots is R = 30. 



Our numerical analysis only depends on the choice of the single parameter n, i.e. the winding 
number. Figure [11] presents the dependence of v and the corresponding energy density S, 
as a function of f, for various values of n. We observe that as the value of n increases, the 
width of the vortex core increases and the energy density radially spreads out, giving the 
characteristic volcano shape. 

5.3 Global Monopoles 

We complete our study of the topological defects that may form in the 2HDM due to the 
spontaneous breaking of the 6 accidental symmetries with the global monopole. This topo- 
logical solution arises from the symmetry breaking of the S0(3)hf symmetry to its subgroup 
S0(2)hf- In spite of being intrinsically unstable, global monopoles may have important cos- 
mological implications, as they can provide a mechanism for structure formation within the 
Universe [50j. 

As with our previous topological defect studies, we do not analyse the cosmological 
implications of the 2HDM's global monopole, but focus on presenting an overview of possible 
solutions. 

5.3.1 SO(3)hf Global Monopoles 

Our analysis in Section l33] has shown that a S0(3)HF-invariant 2HDM potential can exhibit a 
vacuum manifold containing non-contractible 2-spheres, provided a neutral vacuum solution 
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Figure 12: Plots of v{r) and the dimensionless energy density £{r) for the S0(3)hf invariant 
potential. The cut off radius for this plot is R = 20. 



for a global minimum exists, such that (fi)^ + (f^)^ = v^-^. This scenario is only possible in 
the S0(3)HF-invariant 2HDM, for a singular matrix N^,^ [cf. Section [3.3.1) . 

We may seek a time-independent spherically symmetric global monopole solution for 
the S0(3)hf symmetry, by making use of the following ansatz for the two Higgs doublets: 

Mr^x) = -7k\ , , I ' 02(r,X,V') = 4^ ( , , ., I ' (5-35) 

yz I v{r)sinx I \ t'(r)e*''' cos x 

where the coordinate r describes the space radially outward from axis of symmetry of the 
monopole, x is an azimuthal angle and ip is a. polar angle. Using this ansatz, the energy per 
monopole is 

POO 

E = An r'^dr £(0i,02) , (5.36) 
Jo 

where the energy density for the S0(3)hf invariant 2HDM is given by 

Sir) = ^ (^) ' + ^2<rr - l^lvirf + ^XMr)' + Vo , (5.37) 

For this type of energy density, the integral of fl5.36p is linearly divergent in r, so we truncate 
the region of integration from [0, oo) to [0,-R], where is a cut off radius |4lj. Our study 
gets considerably simplified, if we rescale f l5.36p to become dimensionless and introduce the 
dimensionless quantities 

f = LLir , v(r) = , E = . 5.38 
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We then define rj = and require that r] = vsu- Under these rescahngs, the boundary 
conditions on the vacuum field v{r), which follow from (13.371) . become 

lim'O(f) = , lim v{f) = 1 . (5.39) 

f— >0 f—^oo 

These conditions force v{f) to be regular for all values of f and ensures the dimensionful field 
f (r) approaches its VEV in the limit r — )■ oo. The equation of motion for the rescaled field 
v{f) is 

(fv ^ 2dv _ ^/ X + ^ + (5 40) 

As with the majority of monopole studies, we rely on gradient fiow techniques to numerically 
solve the above differential equation. In Figure [121 we present the single solution for an 
S0(3)hf global monopole, by displaying the f-dependence of the vacuum field v{f) and its 
respective dimensionless energy density E. 



6 The U(1)y- Violating 2HDM 

In this section we discuss the application of our Majorana scalar-field formalism to 2IIDM 
potentials which are not restricted by the U(1)y hypercharge group. Even though such 
potentials may not be viable within the context of the SM, they may still be realized in 
models describing cosmological infiation [HI US]. Furthermore, we classify all possible 15 
symmetries that may occur in a general U(l)Y-violating 2IIDM potential, within the 6- 
dimensional bilinear field space. 

If conservation under some U(1)y hypercharge group is lifted from the theory, then 
additional SU(2) gauge-invariant bilinears can, in principle, be present in the 2HDM poten- 
tial, such as (j)Jia'^(j)2 and its Hermitian conjugate, —02^(7^0^ Counting the number of real 
independent parameters, the resulting potential would have 6 bilinear mass terms and 20 
quartic terms. Its explicit analytic form is given by 

V= - /i2(0202) - ml^{4>\(p2) - mt2(020i) - ml^{4>li(^'^<p2) + m;^(0^iaVD 

+Ai(</.l0i)2 + A2(402)' + A3(0l0i)(402) + \M\<t^2M^i) + y (0102)' + f (40i)' 

+ A6(0l0l)(0l02) + A*(0l0i)(0^0i) + A7(0502)(0l02) + \*M<P2M<P1) 

+A8(0l0i)(0^*a202) - A*(0l0i)(05*aVn + \M'p2){'PW'p2) - A;(0502)(05^orVt) 

+ Alo(0l02)(0r^fT202)-Alo(050l)(05«fT20*) + An(050l)(0V02)-Ati(0l02)(05^a20*) 

+^{f.^<r'ct.2f + ^{<t^\^<y^lf. (6.1) 
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We note that the quartic couphngs Ai 2,3,4 are real and A5,6,...,i2 are complex. 

In order to account for the additional bilinear and quartic terms that occur in the 
U(1)y- violating 2HDM potential, we need to promote the 4-vector R'^ in f l2.15p into the 
6-vector R^, with A = 0, 1, 2, 3, 4, 5. The individual components of R^ read: 



R^ 



\ 



^l02 



401 



402 



\ 



— z 



(6.2) 



As with the 4-vector R'^, we can construct the 6-vector R^ using the 8-dimensional complex 
multiplet $ as R^ = To determine the structure of S^, we start again with the 

general GL(8, C) covariant (and SU(2)L-invariant) ansatz 



2^ ap^ 



a'' 



(6.3) 



The particular form of S^^^ is now only constrained by the Majorana condition on S^, 



namely (E 



C "^E^C, in close analogy with (IB.13p . In terms of the tensor S^^o, the 



Majorana condition requires that 



al3 



(6.4) 



Only 6 elements of S^^^^ survive this constraint: S^qq, ^ q^^, ^ Q3, ^ ^ 22 '^^^"^ ^ 32- 



„„, ^^01' ^^03' ^^125 ^^22 ^^'^ ^^32- Hence, 
the six components of the 6-vector compatible with the Majorana condition have the 
tensorial structure 



cr'- 



O2 



Oo 



a 



O2 la' 
^ O2 



a 



O2 



-cr 



—a 
O2 



(6.5) 



Comparing (16. 5p with (IB.20p . we notice that the imposition of the U(1)y hypercharge sym- 
metry on the SU(2)-invariant potential restricts E^'^ = Og and so effectively reduces R^ to 
R'^, as it should. 
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In the absence of the U(1)y hypercharge symmetry, the transformation matrix M no 
longer sphts into two distinct parts, but takes on the general form as determined in fl2.25p . 
Under a SU(2)L-invariant reparameterization-group transformation M G GL(4, C) of the 
scalar-field multiplet $, with M* = CMC [cf. 02.241) ] . the 6- vector R"^ transforms as 

rA ^ R'A _ e'^/^^eR^ , (6.6) 

where e"^ = det[M'''M] > is a real scale factor and A\ is related to the transformation 
matrix M by 

e"/^A\S^ = M^E^M. (6.7) 

Note that the matrix A^g is an element of S0(l,5). This last fact may be verified by defining 
E = (S°, — S^'^'^''^'^), in direct analogy with a'^ = (cr°, — cr^'^'^), and checking the Clifford 
algebra: 

SAs^ + S^S^ = ir/^^Is, (6.8) 

where Ig is the 8-dimensional identity matrix and r]^^ = diag (1,-1,— 1,-1,— 1,-1) is the 
respective metric for the (1 + 5)-dimensional Minkowski fiat space. As a byproduct of (16. Sp . 
we obtain that 

trfS^s''] = 2r/A^. (6.9) 
The latter can be used to compute A^g as 



A^B = ^e-'^/^r/Bctr[MtS^MS^] . (6.10) 

With the aid of the newly introduced 6-vector R^, the potential of (16. ip can be written 
down in a quadratic form similar to (12. 5p : 

V = -1 MaR^ + \ LabR^R"" , (6.11) 

where the 6-vector Ma containing the mass terms and the 6x6 quartic coupling matrix Lab 
read: 
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MA = (/if+At|, 2Re{mf2), -2Im(m^2) , ^t^ - ^2 ^ 2Re{ml^) 

( Ai + A2 + A3 Ro(A6 + A7) 

Re(A6+A7) A4 + Ro(A5) 
~Im(Ae+A7) -Im(A5) 

Ai — A2 Ro(A6 — A7) 



-2Im(m,34) 



(6.12a) 



Lab 



Re{A8 + Ag) 
V -Im(A8 + Ag) 



Re(Aio + An) 
-Im(Aio + All) 



-Im(A6 + A7) 

-Im(A5) 
A4 - Rc(A5) 
-Im(A6 - A7) 
-Im(Aio - All) 
~Rc(Aio - All) 



Ai — A2 
Re(A6 - A7) 
-Im(A6 - A7) 
Ai + A2 — A3 
Rc(As — Ag) 
-Im{A8 — Ag) 



Re(A8 + Ag) 
Re(Aio + All) 
-Im(Aio - All) 
Ro(A8 - Ag) 
Rc(Ai2) 
-Im(Ai2) 



-Im(A8 + Ag) \ 
-Im(Aio + All) 
-Rc(Aio - All) 
-Im(A8 - Ag) 
-Im(Ai2) 
-Re(Ai2) 



{6.12b) 



Note that in the U(l)Y-symmetric hmit, Ma — J- and Lab — ^ L^^/, whereas the elements 
of Ma and Lab vanish for the components A, B = 4, 5. 

We may now use an approach analogous to [22|, in order to identify all possible acci- 
dental symmetries that could take place within a general U(l)Y-violating 2HDM. Requiring 
that the kinetic terms remain invariant under GL(8, M) scalar-field transformations, we are 
restricted to consider unitary rotations U G U(4) in the $-space, subject into the Majorana 
constraint: U* = CUC. These Majorana-constrained U(4) transformations induce orthog- 
onal rotations S0(5) C S0(1, 5), which act on the 'spatial' components A = 1, 2, ... 5 of the 
6-vector R^. In detail, we may classify all possible symmetries derived from S0(5), which 
include SO (5) and its proper, improper and semi-simple subgroups. If Z2 is the reflection 
group for one of the spatial components of R"^, we may now list all the symmetries starting 
from the larger and going to the smaller group. In this way, the symmetries may be grouped 
into the following five categories: 



L S0(5); 

II. 0(4) ®Z2; S0(4); 

III. 0(3)® 0(2); SO(3)®(Z2)'; 0(3) ® Z2; S0(3); 



(6.13) 



IV. 0(2)®0(2)®Z2; 0(2)® 0(2); 0(2) ® (Z2)^ S0(2) 
0(2) ®Z2; S0(2); 



fZ2)2; 



V. 



[Z2)^; (Z2)^ 



Note that all the symmetry transformations have determinant equal to +1. With this re- 
striction, we get 15 distinct symmetries that could act on a general tree-level U(1)y- violating 
2HDM potential. Moreover, the above classification in (16.131) contains the U(1)y group. More 
explicitly, the six accidental symmetries reported in the literature are: the first symmetry un- 
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Symmetry 




m^2 


77134 Ai 


A2 A3 


A4 Re A5 


Ag — A7 


Ag — Ag 


Z2 


- 








- 


- 








U(1)PQ 



















S0(3)hf 


- ^^i 







Ai - 


2Ai - A3 








CPl 




Real 


Real - 






Real 


Real 


CP2 


- i^i 





Real - 


Ai - 







Real 


CPS 


- i^i 





Real - 


Ai - 


2Ai — A34 









Table 12: Parameter relations in the general U {l)Y-violating 2HDM potential that result from 
the imposition of the six accidental symmetries, in the diagonally reduced basis ImAs = 0, 
Aio = All = and ImAi2 = [cf. fl6.14p /. The quartic coupling ReAi2 remains uncon- 
strained by the six considered HF/CP symmetries. Finally, a dash indicates the absence of a 
constraint. 



der Category III and the first 5 symmetries under Category IV, i.e. 0(3)00(2); 0(2) 00(2) (g) 
Z2; 0(2)®0(2); 0(2)® (Za)^ SO(2)®(Z2)^ 0(2)®Z2. In Table [IH we show the parameter 
restrictions of these six HF/CP symmetries for the full U(l)Y-violating 2HDM potential, as 
these are realized in a specific basis where the spatial part of Lab (with A, B = 1, 2 ... 5) is 
made diagonal by an SO (5) rotation. In such a diagonally reduced basis, we have 

ImAs = 0, Ae = A7 , Ag = Ag , Aio = An = , ImAi2 = . (6.14) 



Given the classification in (16.131) . we observe that symmetries higher than 0(3), which 
contain the U(1)y group, can still occur. For instance, one such symmetry is S0(5), which 
is obtained when 2Ai = 2A2 = A3, fif = and all other parameters vanish. The symmetry 
SO (5) is equivalent to 0(8) [21] in the real field space and includes the gauge-group rotation 
SU(2)l®U(1)y. In the extended bilinear R^-space, S0(5) breaks down to S0(4) or 0(4) x Z2, 
giving rise to four pseudo-Goldstone bosons, as it should be. Notice that within the SU(2)l 
and U(1)y constrained bilinear formalism, it is not possible to clearly make the distinction 
between the S0(3)jjp symmetry and the possible higher HF/CP symmetry S0(5). 

Another interesting example is the symmetry S0(4), which is obtained from a U(1)y- 
and Za-invariant 2HDM potential, with the additional constraint that A4 = A5 = 0. This 
model is equivalent to the model 0(4) ® 0(4) [31J in the scalar-field space, where the second 
0(4) describes the gauge group rotations. The symmetry SO (4) breaks into SO (3), giving 
rise to three pseudo-Goldstone bosons. Again, this breaking scenario cannot be distinguished 
within a SU(2)l and U(1)y constrained bilinear formalism, and can be easily confused with 
the CP3 symmetry. In Table [T31 we display the 7 additional accidental symmetries that may 
occur in a U(l)Y-invariant 2HDM potential, along with parameter restrictions obtained in 
the diagonally reduced basis [cf. (16.141) ]. Note that all symmetries lead to CP-invariant scalar 
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Symmetry 




771^2 Al 


A2 


A3 


A4 


Re A5 


Ae — A7 


S0(5) 


- f^i 





Al 


2Ai 











0(4) X Z2 


- 





Al 


- 











S0(4) 




















0(3) X 0(2) 


- f^i 





Al 


2Ai 










S0(3) X (Z2)" 


- f^i 





Al 






A4 





0(3) X Z2 


- f^i 


Real - 


Al 






A4 


Real 


S0(3) 




Real - 








A4 


Real 



Table 13: Parameter relations in the general 1]{1)y -invariant 2HDM potential that result 
from the imposition of the additional accidental symmetries shown in Categories I, II and 
III of (16.131) . in the reduced basis ImAs = and Xq = A7 [cf. (16.141) /. A dash indicates the 
absence of a constraint. 

potentials. Further details of these additional HF/CP symmetries will be given elsewhere. 

7 Conclusions 

Unlike the SM, the 2HDM has a rich landscape of discrete and global symmetries, whose spon- 
taneous breaking may lead to non-trivial topological solutions. In this paper, we have taken 
the first step towards analyzing a number of generic symmetries for their resulting vacuum 
topology within the 2HDM. For definiteness, we have considered the three HF symmetries: 
Z2, U(l)pQ and S0(3)hf, and the three CP symmetries: CPl, CP2 and CP3 (cf. Table 2). 
In order to study the vacuum topology of these six symmetries, we have introduced a Ma- 
jorana scalar-field formalism based on two subgroups of GL(8,C), where the HF and CP 
transformations may act on a single scalar-field multiplet representation. 

Using Sylvester's criterion, we have derived the general conditions in order to have a 
convex, stable and bounded-from-below 2HDM potential. Given these convexity and stability 
constraints, we have solved analytically the minimization conditions of the scalar potential, 
by making use of the Lagrange multiplier method. We have thus obtained all two non-zero 
solutions for the neutral vacuum expectation values of the Higgs doublets for the aforemen- 
tioned six HF and CP symmetries, in terms of the gauge-invariant parameters of the theory. 

In order to identify the nature of the topological defects associated with the sponta- 
neous symmetry breaking for each of the above six symmetries, we have studied the homotopy 
groups of the resulting vacuum manifold after spontaneous symmetry breaking. In partic- 
ular, we have found the existence of domain walls from the breaking of Z2, CPl and CP2 
discrete symmetries, vortices in models with broken U(l)pQ and CP3 symmetries and a global 
monopole in a model with SO(3)jjp-broken symmetry. We have then studied the topological 
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defect solutions numerically, as functions of the potential parameters of the 2HDM. We have 
given numerical examples for each topological defect, showing the energy of the defect for 
typical situations. 

As we have explicitly demonstrated in Section El our Majorana scalar-field formalism 
can be applied to identify 7 further accidental symmetries in the 2HDM potential, which 
include the maximal symmetries 0(8) and 0(4) (8> 0(4) in the real field space [31]. These 
symmetries remain undetected by the constrained SU(2) bilinear field approach considered 
so far in the literature. 

Our Majorana scalar-field formalism can also be used to study more general scalar 
potentials which are not constrained by the U(1)y hypercharge symmetry and can realize 
a maximal number of 15 distinct symmetries. Such 2HDM potentials may not be directly 
related to the observable SM gauge group, but may form an independent hidden sector, 
as it is, for example, the case in supersymmetric theories of hybrid inflation [HI US]. The 
formation of topological defects, such as domain walls, cosmic strings, monopoles, or textures, 
through the spontaneous symmetry breaking of global, local or discrete symmetries may 
have important implications for the analysis of the cosmological data. It would be therefore 
interesting to analyze the cosmological constraints on the fundamental parameters of the 
2HDM, using the formalism and the computational framework developed in this paper. 
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Appendices 



A (j^ Matrix Identities 

Here we list a number of useful identities for the matrices a'^ = (a", cr^'^'^), where o"° = I2 
and 0"^'^'^ are the standard Pauli matrices. These identities are used in Appendix [B] to 
derive the explicit form of S^. Under transposition and complex conjugation, the individual 
components of a'^ transform as 





)^ 






)* 




a': 








r 






)^ 






)* 


= -a^ 




)^ 






r 


= a\ 



Hence, the above identities may be cast into the more compact form: 

i^'f = , (A.la) 

{an* = {6^rX , (A.lb) 

with 

^ diag(l,l,±l,l). (A.2) 
We will also frequently use the sandwich products 

aV^a^ = {Ji)y , (A.3a) 

aW = (J2)>% (A.3b) 

aV^a^ = (3,)^^ , (A.3c) 

where the tensors Ji,2,3 are defined as 

(Ji)t = diag(l, 1,-1,-1) , (A.4a) 

(J2)t = diag(l, -1,1,-1), (A.4b) 

(J3)t = diag(l, -1,-1,1). (A.4c) 

Finally, it is interesting to note the identity 

ihnS-)^ = < • (A.5) 
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B The Form of and the Transformation Matrices 



In order to derive the explicit form of T,^ in GL(8, C), we start with the following general 
ansatz: 

E^^ = S%cT°®a^ (B.l) 

where we have suppressed the SU(2)l gauge-group space for convenience. Then, we need to 
apply two constraints to determine the tensor coefficients S'^^: the U(1)y constraint and the 
Majorana constraint. 

B.l The U(1)y Constraint on E^^ 

Under a U(1)y transformation, the 4-component multiplet $ defined in fl2.13p transforms as 
follows: 

$' = Uy$ , (B.2) 

where 

Uy = e'^'^-'^-"^ = diag (e^^^ e^', e'^^^ e'^^^) = B^a' ® a° , (B.3) 

with 

= [cos (Y^), 0, 0, isin(Y^)] . (B.4) 

Invariance of the 4-vector = [cf. fl2.15p ] under a U(1)y transformation implies the 

following double equality constraint on S'^: 

S'' = UyS^Uv = UyS'^Uy . (B.5) 

Given the ansatz of in (IB.ip . the above double constraint gets translated into: 



E'^ = U:^S"Uy = S%B:Ba [{aya'^cr^] ® , (B.6a) 
S'^ = UyS^U^ = S'^^B^B* [(tV°((t^)*] ® . (B.6b) 



Using the identity (lA.lbp . the above two relations can be rewritten as 



U^S^Uy = S%B:B4(5_);aV°a^] ®a^ (B.7a) 
UyS^U:^ = S%B,B*[(5_)>V"a^] (B.7b) 
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Substituting the explicit forms of and f lB.7aP and f lB.7bp become respectively: 



E'^ = S% (^cos2(Y0)(T" + sin2(Y0)(J3)"^a^ + isin(Y^)cos(Y0) [fT",^^]) ® , (B.8a) 
W = S';^ (^cos2(Y^)a" + sin2(Y^)(J3)"^a^ -isin(Y^)cos(Y^) ^ . (B.8b) 

Evidently, in order that the above two constraints are satisfied, the commutator term must 
vanish, i.e. 

[cr°,cT3] = . (B.9) 
This can only happen for the choices a = 0, 3, implying that 

^'ip = = , (B.IO) 

independently of the Lorentz indices /i and f3. As a consequence, the U(1)y constraint leads 
to the block diagonal form for the matrix T,^: 

B.2 The Majorana Constraint on 

The Majorana condition fl2.18p on the scalar multiplet $ gives rise to another important 
constraint on the form of E'*. Specifically, the condition fl2.18p implies the invariance of 
vector defined in f l2.15p under charge conjugation. Thus, when $ — C $*, R'^ transforms 
as 

RA* = ^^W^ ^ R^ = $TctE^C$* = <l>"fC^(E^)^C*$ . (B.12) 
Requiring that R'^ = R^ yields the Majorana constraint: 

(E'^)^ = C^^E^C. (B.13) 
For the general ansatz (IB. 10 . the last constraint is equivalent to 

K^i^'^V ® i^^f = S% (aW^) ® . (B.14) 
Employing the identities of Appendix [XJ we obtain the constraining equation on E'^^: 

= ^'x,viiS.r, . (B.15) 
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Assuming that S'^ has the U(l)Y-invariant form f IB.lip and using the identity f lA.Sp allows 
us to express S'^^ as follows: 



Sy5_r^, fora = 
-^%pi^-Yp . for a = 3 



(B.16) 



From this last expression, we find that the two non-zero parts of the E'^^ tensor are then, in 
general, proportional to the following matrices: 

oc {6^r^ + {6_r^, (B.17a) 
S'^sp « {6^r,-{6-r,. (B.lTb) 

This can be written down in the covariant form: 

E% = 8.^is+r^+h^is.r^, (B.18) 

where the vectors a.a and ha are defined as 

a„ = J (1,0, 0,-1) , (B.19a) 
K = ^(1,0,0,1) . (B.19b) 

Implementing all the above results, the U(l)Y-invariant vector compatible with the Ma- 
jorana constraint takes on the simple form: 

1 f a^' 

E'^ = - I . (B.20) 

2 (a'^V ' 



B.3 The Majorana Constraint on GL(8,C) 

It is interesting to discuss the reduction of the GL(8, C) group under the Majorana constraint 
M* = C M C for HF symmetries, where M = M^,,a^ (g) a" (with M^^, G C) becomes a 
general member of GL(4, C) after suppressing the SU(2)l gauge group space. The Majorana 
reduction, M, pertinent to CP transformations is analogous and will not be repeated here. 
Applying the Majorana constraint on M, we obtain the expression 

M* = M;, {a^^y ® {an* = M^, [a^ ® a") (a^ ® a") [a' ® a°) . (B.21) 
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We may now use the so-called mixed-product identity: (A (g) B) (C ® D) = (AC) (g) (BD) and 
the identity (]A.lb|) . in order to rewrite (lB.2ip as follows: 



a 



(B.22) 



Further use of the sandwich products given in Appendix |A] implies 
which translates into the constraining equation: 

Solving this last equation term by term results in the following constraints: 



(B.23) 



(B.24) 



Moo 


= Mo*o 


Moi 


= Mo*i 


M02 


= -Mo*2 


Mo3 = 


Mo*3 


Mio 


= -Mi*o 


Mn 


= -M^i 


M12 


= Mi*2 


Mi3 = 


-Mr3 


M20 


— —M* 


M21 


= -M2*i 


M22 


= M*2 


M23 = 


-M2*3 


M30 


— —M* 


M31 


= -M3*i 


M32 


= M|2 


M33 = 


-M3*3 



Hence, from the 32 independent parameters of M, half are eliminated by the Majorana 
condition. The resulting 16 free parameters generate a group which is isomorphic to GL(4, R) 
acting on a complex four-dimensional vector space. 



C Trace and Determinant Relations for N^j, and L^^^ 

Relations involving the traces and determinants of N^,^ and L^^, play an important role in the 
calculation of the VEVs of the Higgs doublets and in the derivation of stability and convexity 
conditions for the 2HDM potential. 

To facilitate our presentation, we use the shorthand notation N = N^j,, L = and 
7] = Tj^y to represent the 2-rank tensors as 4 x 4 matrices. We also assume the standard 
multiplication law between matrices, e.g. (N^)^,^ = N^q,Nq,^, (Lr^)^,^ = h^aVau etc. In the 
above notation, the determinant of N may be written as 

det[N] = det[L-Cr?] , (C.l) 

which can be calculated by the following determinant-trace identity: 

det [N] = — {tr^ [N] - etr^ [N] tr [N^] + Str^ [N^] + 8tr [N] tr [N^] - 6tr [N^] } . (C.2) 
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The trace relations between N and L are found to be 



tr [N] = tr [L] + 2C , (C.3a) 

tr [N^] = tr [L^] - 2Ctr [Lt^] + 4C^ , (C.3b) 

tr [N^] = tr [L^] - 3Ctr [h'^rj] + SC^tr [L] + 2C^ , (C.3c) 

tr [N^] = tr [L^] - 4Ctr [L^r)] + S^tr [L^] + 2C^tr [L77L77] - 4C^tr [Lrj] + 4C^ . (C.Sd) 

Thus, the determinant of N is given by 

det[N] = -(^-AC'-BC-CC-D, (C.4) 

where 

^=-tr[L77], (C.5a) 

B =tT[L^]-hT%r]]+2tT[L^r)] + ^tr[Lr]Lr]]-tT[L] (2tr[Lr7] +tr[L]) , (C.5b) 

C=-tr[L3^]+tr[L] (trfL^Tyl+trfL^D + ^tr^ (tr[LVtr'[L]) -^tr^M -^tr[L^] , (C.5c) 

D=-det[L]. (C.Sd) 

Notice that the coefficients A, C and D are entirely expressed in terms of traces of powers 

of L and the determinant of L. These latter expressions depend explicitly on the quartic 
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couplings of the 2HDM potential as follows: 



tr[L] 
tr [L^] 



2Ai + 2X2 + 2A4 , 

AXl + AXl + 2X1 + 2A4 + 8IA5I' + 4|A6|' + 4|A^ 



(C.6a) 
(C.6b) 



tr [L^] = 8X1 + 8^2 + 6(Ai + X2)Xl + 2X1 + 24A4IA5I' + 12Ai|A6|' + I2A2IA7P 
+ 12A3(i?6^7 + W + 6(A4 + 2R5){Rl + i??) + 6(A4 - 2i?5)(/| + I7) 



tr [L^] = 16Xt + l6Xt + 2A3 + IGiXl + Xl + AiA2)A^ + 2A^ + WlX^W^Xl + 2\X5\^) 
+ 8|A6|'(4A? + 2A1A4 + A^ + A^ + 4IA5I' + lAel') 
+ 8|A7p(4A^ + 2A2A4 + A^ + A^ + 4IA5P + lAyp) 

+ 32(Ai + A4) [R^iRl - II) + 2hRd^] + 32(A2 + A4) [R^{R^j - I^) + 2hR7h] 
+ 16(i?6i?7 + hhf + 16A3(2Ai + 2A2 + X^){ReR7 + hh) 

+ 32A3 [R^iRf^Rj - hh) + h{Rdi + hRi)] , (C.6d) 
tr [1.7]] = 2X3 - 2A4 , (C.6e) 
tr [L^T]] = 4(Ai + A2)A3 - 2A2 - 8IA5I' - 2{Re - R^f - 2{h - h? , (C.6f) 
tr [L^ = A{2Xl + 2X1 + 2A1A2 + lAgp + lArHAg - 4Ai|A6p - 4A2IA7P 

+ 4(2Ai + 2A2 - X3){ReR7 + hh) + 2A^ - 2A^ - 24A4IA5I' 

- 4(A4 + 2R,){Rl + i?2 - R,R,) - 4(A4 - 2R,){ll + - hh) 

+ 8hihR7 + Reh - 2R7h - 2Reh) , (C.6g) 
tr [Lr/Lr/] = 8A1A2 + 2X1 + 2A4 + 8IA5I' - 8(i?6^7 + hh) , (C.6I1) 
det[L] = (4A1A2 - XDiXl - 4IA5I') - 4A4(Ai|A7r + A2|A6r) + A\Xe\'\X7\' 

- A{ReR7 + hh)^ + 8A1 [R,iR^, - /2) + 2hR7h] 

+ 8A2 [R^iRl - II) + 2hRd&] + 4A3A4(/?6i?7 + hh) 

- 8A3 [R^{RqR7 - hh) + HhRi + ^6/7)] ■ (C.6i) 



To find the values for the Lagrange multiplier C, that lead to a singular N matrix with 
det[N] = 0, we need to solve a quartic equation. To do so, we first apply the standard linear 
transformation to C, 



which enables one to reduce the quartic order polynomial of f lC.4p to the incomplete quartic 



+ 24h{R6h + R7h) 



(C.6c) 




(C.7) 



equation 



+ ap^ + /3p + 7 = , 



(C.8) 
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where 





B- 


3^2 


a = 


8 


P = 




AB 


8 


2 




D - 


AC 


7 = 


4 



c 

/l2 ] 



16 256 

In terms of the quartic couphngs Ai,2,...7, the coefficients a, f3 and 7 are given by 

a = -4AiA2-i(A3 + A4)'-4|A5|' + 4(i?6i?7 + /6/7) , 

(3 = (4|A5p-4AiA2)(A3 + A4)+4Ai|A7|' + 4A2|A6r 

— SR^^RqRy — IqIj) — SI^IIqRj + -Rg-^t) , 

7 = 16AiA2|A5p + 4(^3 + A4r-(A3 + A4)'(AiA2 + |A5p + i?6i?7 + /6/7) 
16 

+ 2(A3 + A4)(Ai|A7p + AalAeH - 4|A6nA7|' + AiR^R^ + hhf 

- 8A1 [R^{R\ - J2) + 2hRjh] - 8A2 [R^iRl - II) + 2hRd^] 

+ 4(A3 + A4) [R^{RgR7 - hh) + HRdi + /6i?7)] • 



(C.9a) 
(C.9b) 
(C.9c) 



(C.lOa) 
(C.lOb) 



(C.lOc) 



The analytical solutions to the incomplete quartic equation can now be found by making 
using of the Descartes-Euler method. To this end, we first construct the cubic resolvent 
equation of flC.8p . which is 



+ 2ax^ + (a^ - 47)0; - /J^ 
whose roots are determined by the standard formulae: 

„2 







(C.ll) 



Xi 



X2 



X3 



6 2a' + 247 _ 2a 
6 ^ 35 3" 



{l + iVS)5 {1 - iVS){a' + 12-f) 



12 



36 



{l-iV3)5 (l + iy3)(a2^127) 



12 



36 



2a 

T 

2a 



(C.12a) 
(C.12b) 
(C.12c) 



where 



(J3 = 80^-28807+108/32 



+12y/-A8a^-f + 384^2^2 _ 753^3 + Ua^/S'^ - 4320/327 + 81/3^ . (C.13) 
Having thus obtained the cubic roots xi,2,3, the four roots Ci, 2,3,4 of the original quartic 
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equation det[N] = are then given by 



/■ _ 

u — 


1 

~2 


,V^l + V 




A 
4 ' 


(^0. i4a j 


C2 = 


1 

-2 




■ y^) - 


4 ' 


(C.14b) 


C3 = 


1 

-2 


;-v^ + 




~ 4 ' 


(C.14c) 


C4 = 


1 

-2 




+ V^) 


~ T ■ 


(C.14d) 



D Inverting the Transformation Matrix Relations 

It would be useful to give the relations between the transformation matrices M and the 
S0(l,3) matrices A'^, by assuming that the scale factor is e*^ = 1. 

As was discussed in Section 12.2^ the general matrix M may describe both the HF and 
CP transformations by the matrices M+ and M_, respectively, which contain the reduced 
two-by-two matrices T-t G SL(2,C). Following [75], we first note that 

o^T^cy^ = 2tr[T±]or°, (D.l) 

where = (cr°, — cr^'^''^). On the other hand, contracting f l2.3ip and (12.351) from the RHS by 
S'^ = diag [a'^ , (o''^)"'"] yields the relations 

(A±)t (^±) A^'^M = ^icr.T^^^ . (D.2) 
Making use now of the identity (ID.ip . we can solve for Tj_: 



Tl = {S±)\{A±n cr% . (D.3) 

To remove the a'^- dependence from the RHS of the above equation, we use a relationship 
derived by taking the determinant on both sides of (ID. 21) : 

det[(A±)t V^T^l = det[T^(2tr[T±])] = (2tr[T±])2 det[T^] . (D.4) 
Since det [T±] = 1 for T± G SL(2, C), one arrives at 

2tr[T±] = {det[(A±)t (5±)>V^]}^ . (D.5) 
Here, we have omitted the negative solution from the square root as this is accounted for by 
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the U(1)y invariance of the theory. Thus, one ends up with the expression 

Ti = ^ r A(A±)t ^'^M • (D.6) 

{det [(A±)t 

The determinant in the denominator of the above equation can be calculated using the 
relation: 2det[G] = tr[G]^ — tr[G^], which results in 

= 4 + tr[A±5±]2 - tr[(A±5±)2] - ^e^'^''"(A±)^,(A±)„^(5±)^,(5±)^, . (D.7) 

Here we use the convention e°^^'^ = +1 for the Levi-Civita tensor. We can now use the 
identity 

cr'cr, = + v,'^' - r/^V, + ^e'\^a^ , (D.8) 

to write down the numerator of (]D.6P in the form 

(5±)^,(A±)t = tr[A±<5±]a° + {{6^yi{A^r^ - (5±);(A±)/ + ^e'\,{6^n{A^r^} . 

(D.9) 

Using the representation T-t = (T±)^cr^, the individual components of (T-j.)^ derived from (]D.6P 
are given by 

(T±)o = T^tr[A±5±] , (D.lOa) 
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